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ABSTRACT 


A  new  slender-ship  theory  of  wave  resistance  is  presented.  Specifically, 
a  sequence  of  explicit  slender-ship  wave-resistance  approximations  is  obtained. 
These  approximations  are  associated  with  successive  approximations  in  a  slender- 
ship  iterative  procedure  for  solving  a  new  (nonlinear  integro-dif ferent ial) 
equation  for  the  velocity  potential  of  the  flow  caused  by  the  ship.  The  zeroth-, 
first-,  and  second-order  slender-ship  approximations  are  given  explicitly  and 
examined  in  some  detail. 


The  zeroth-order  slender-ship  wave-resistance  approximation,  r 


(0) 


obtained  by  simply  taking  the  (disturbance)  potential,  4> ,  as  the  trivial  zeroth- 

order  slender-ship  approximation  =  0  in  the  expression  for  the  Kochin  free- 

wave  amplitude  function.  The  classical  wave-resistance  formulas  of  Michell  (1898) 

and  Hogner  (1932)  correspond  to  particular  cases  of  this  simple  approximation. 

The  low-speed  wave-resistance  formulas  proposed  by  Guevel  (1974),  Baba  (1976), 

Maruo  (1977)  and  Kayo  (1978)  are  essentially  equivalent  (for  most  practical  purposes) 

to  the  first-order  slender-ship  low-Froude-number  approximation,  r»^,  which  is  a 

(lj 

particular  case  of  the  first-order  slender-ship  approximation  r  .  Specifically, 

the  first-order  slender-ship  wave-resistance  approximation  r^  is  obtained  by 

approximating  the  potential  4>  in  the  expression  for  the  Kochin  function  by  the  first- 

order  slender-ship  potential  <t>^\  whereas  the  low-Froude-number  approximation  r p 

(1)  (1) 

is  associated  with  the  zero-Froude-number  limit  of  the  potential  <(> 

u  (1) 

A  major  difference  between  the  first-order  slender-ship  potential  <f>  and 


(1) 


resides  in  the  waves  that  are  included  in  the 


its  zero-Froude-number  limit  <p 

(1)  u 

potential  <J>  but  are  ignored  in  the  zero-Froude-number  potential  4> 


(1) 


Results 


of  calculations,  due  to  C.Y.  Chen,  for  the  Wigley  ship  form  show  that  the  waves  in 
the  potential  $ ^  have  a  remarkable  effect  upon  the  wave  resistance;  in  particular, 
they  cause  a  large  phase  shift  of  the  wave-resistance  curve  towards  higher  values 
of  the  Froude  number.  Comparison  of  the  first-order  slender-ship  wave-resistance 
approximation  with  experimental  results  shows  fairly  good  agreement,  in  considerable 
improvement  with  respect  to  the  Guevel-Baba-Maruo-Kayo  low-speed  theory. 


—"•fir 
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1.  Introduction 

The  problem  of  predicting  the  wave  resistance  of  a  ship  in  steady  rectilinear 
motion  in  a  quiescent  sea  is  one  of  the  classical  and  basic  problems  in  ship  hydro¬ 
dynamics.  It  is  also  one  of  the  most  difficult.  As  a  matter  of  fact,  in  spite  of 
considerable  research  activities,  notably  in  the  past  decade,  since  the  famous 
pioneering  study  of  Michell  (1898),  discrepancies  between  values  of  the  wave- 
resistance  coefficient  determined  experimentally  and  predicted  by  various  analytical 
theories  and  numerical  methods  often  are  quite  significant  and  may  even  be  extremely 
large  in  some  cases,  as  is  clearly  apparent  from  the  Proceedings  of  the  recent  Work¬ 
shop  on  Ship  Wave-Resistance  Computations  (1979). 

In  this  study,  a  new  slender-ship  theory  of  wave-resistance  is  presented. 
Specifically,  a  sequence  of  explicit  slender-ship  wave-resistance  approximations  is 
obtained.  These  approximations  are  associated  with  successive  approximations  in  a 
slender-ship  iterative  procedure  for  solving  a  new  (nonlinear  integro-dif f erential) 
equation  for  the  velocity  potential  of  the  flow  caused  by  the  ship.  The  zeroth-, 
first-,  and  second-order  slender-ship  approximations  are  given  explicitly  and  ex¬ 
amined  in  some  detail. 

The  zeroth-order  slender-ship  wave-resistance  approximation,  r^\  is  obtained 

by  simply  taking  the  (disturbance)  potential,  <j> ,  as  the  trivial  zeroth-order  slender- 

ship  approximation  4>^^r0  in  the  expression  for  the  Kochin  free-wave  amplitude  function 

The  classical  wave-resistance  formulas  of  Michell  (1898)  and  Hogner  (1932)  correspond 

to  particular  cases  of  the  zeroth-order  slender-ship  approximation  r^\  Specifically, 

the  Hogner  approximation  is  obtained  by  neglecting  the  waterline  integral  (of  order 
3 

b  in  the  beam/length  ratio  b)  in  comparison  to  the  hull  integral  (of  order  b)  in 
the  expression  for  the  zeroth-order  slender-ship  approximation  to  the  Kochin  free- 
wave  amplitude  function.  The  Michell  approximation,  on  the  other  hand,  corresponds 
to  the  consistent  first-order  thin-ship  approximation  to  the  zeroth-order  slender- 
ship  approximation.  A  simple  example  is  used  to  demonstrate  that  neglect  of  the 
waterline  integral  or/and  use  of  the  thin-ship  approximation  may  have  important 
effects. 

The  first-order  slender-ship  wave-resistance  approximation,  r^  ,  is  obtained 
by  approximating  the  potential  <t>  in  the  expression  for  the  Kochin  function  by  the 
first-order  slender-ship  potential  <J>^^.  In  the  thin-ship  limit,  the  potential 
becomes  the  classical  Michell  potential;  more  precisely,  the  Michell  potential  cor¬ 
responds  to  the  first-order  approximation  in  a  thin-ship  expansion  of  the  potential 
♦  In  the  zero-Froude-number  limit,  the  potential  <f>^  becomes  the  potential 


NOTE 


This  report  is  based  upon  three  previous  reports  entitled: 

"Potential  Theory  of  Steady  Motion  of  Ships,  Parts  1  &  2" f  /" 

MIT  Dept,  of  Ocean  Engineering  Rep.  No.  78-4,  Sep.  1978,  43  pp. 

"Part  3:  Wave  Resistance"  MIT  Dept.  Ocean  Eng.  Rep.  No.  78-5,  Nov.  1978,  26  pp. 
"Part  4:  Low-Froude-Nuraber  Approximations"  MIT  OE  Rep.  No.  79-1,  May  1979,  40  pp.* 
A  large  part  of  the  material  included  in  this  report  has  thus  been  presented 
previously.  In  particular,  the  new  integro-dif ferential  equation  (5.1)  for 
determining  the  velocity  potential  of  the  flow  caused  by  the  ship,  upon  which 
the  present  slender-ship  theory  is  based,  was  first  presented  in  a  manuscript 
entitled:  "The  Potential-Flow  Problem  in  the  Theory  of  Steady  Motion  of  a  Ship" 
that  was  submitted  for  publication  to  the  Journal  of  Ship  Research  in  January  1979 
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which  is  a  first-order  slender-ship  approximation  to  the  zero-Froude-number 

(double-hull)  potential  <t>n.  It  is  shown  in  Noblesse  and  Triantafyllou  (1980) 

(1)  U 

that  the  potential  <D  provides  a  fairly  good  approximation  to  the  potential  <j>_ 

0  (n  u 

for  longitudinal  motion  of  a  slender  body.  In  particular,  the  potential  <j>g  is 

proportional  to  the  potential  <J)q  for  ellipsoidal  hull  forms. 

The  wave- resistance  approximation  associated  with  the  use  of  the  zero-Froude- 

number  potential  as  an  approximation  to  the  potential  <j>  in  the  expression  for 

the  Kochin  function  is  referred  to  as  the  low-Froude-number  approximation,  r,  . 

This  wave-resistance  approximation  is  identical  to  the  (closely-related  and  essen- 

tially-equivalent)  low-speed  wave-resistance  formulas  proposed  by  Guevel,  Vaussy, 

and  Kobus  (1974),  Baba  (1976),  Maruo  (1977),  and  Kayo  (1978).  The  wave-resistance 

approximation  associated  with  the  first-order  slender-ship  approximation  (j>g^  t0 

the  zero-Froude-number  potential  <j>.  thus  corresponds  to  a  first-order  slender-ship 
(1)  ° 

approximation,  r^  ,  to  the  Guevel-Baba-Maruo-Kayo  (GBMK)  low-Froude-number  approx¬ 
imation  r^. 

The  low-Froude-number  wave-resistance  approximation  r.  ,  the  first-order 

(1) 

slender-ship  low-Froude-number  approximation  r^  ,  and  the  zeroth-crder  slender-ship 
approximation  r^  (associated  with  the  approximations  4>=<pQt  >  an(*  <+>=$^^  =  0, 

respectively)  are  compared  to  one  another  in  the  particular  case  when  the  "ship" 
hull  is  a  vertical  cylinder  with  elliptic  waterline.  For  this  particular  case,  the 

ClF 


approximations  r  ,  r»  ,  and  rp  ,  actually  are  proportional  to  one  another. 

(CO  '-*2  Ml  2  4 

Specifically,  we  have  rv  =  r^/ (1+b)  ,  and  r^  =  r£F(l+2b)  /(1+b)  ,  where  b  is  the 

beam/length  ratio  of  the  elliptical  cylinder.  The  relative  errors,  e  ^  and  say, 


.(0) 


associated  with  the  approximations  r  and  t'c 


CD  .  .  ,,  ,  (0)  ,  (0)w 

,  and  defined  as  c  =(r^-r  )/r^ 

and  e  ^■(r^-r^^)  /  r£F  >  can  then  readily  be  determined  in  terms  of  the  beam/length 
ratio  b.  It  may  thus  be  found  that  we  have  .174  and  .016  for  b=  .1,  and 

£  ^  .306  and  =  .055  for  b=.2.  The  first-order  slender-ship  low-Froude-number 

approximation  r^^  thus  differs  from  the  low-Froude-number  approximation  r^  by  only 
a  few  percent  for  a  thin  vertical  elliptical  cylinder. 

Calculations  by  Chen  and  Noblesse  (1980)  for  the  case  of  a  vertical  cylinder 
with  a  waterline  in  the  shape  of  an  ogive  have  shown  that  differences  between  the 
approximations  r»^  and  r»  are  somewhat  larger  than  for  an  elliptical  cylinder  with 


same  beam/length  ratio,  but  remain  small,  of  the  order  of  a  few  percent  for  typical 
values  of  b.  This  numerical  study  also  includes  results  of  calculations  for  ellip¬ 
soidal  hull  forms;  these  show  that  the  approximations  and  r ^  are  practically 

indistinguishable  for  values  of  the  beam/length  ratio  equal  to  .1  or  .2.  Finally, 

g  L 

figure  A,  due  to  C.Y.  Chen  ,  for  the  Wigley  ship  form  show  that  the  wave-resistance 
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c  (1 ) 

curve  corresponding  to  the  approximation  r^  is  fairly  close  to  the  points 

determined  by  averaging  the  numerical  results  obtained  by  several  researchers^  on 

the  basis  of  the  GBMK  low-speed  theory. 

The  first-order  slender-ship  low-Froude-number  approximation  r!^ ,  associated 

(1)  . 

with  the  approximation  <}>  =<{> q  ,  is  a  particular  case  of  the  first-order  slender-ship 
approximation  r^,  which  corresponds  to  the  approximation  as  was  noted 


previously.  A  major  difference  between  the  potentials  iji 

i  ii.  i  i — >  .  .  n  /i  1  ■  ■  i — L  .4  i  t.  .  Vi  r — i  n  , .  r\  .  .  1  ^  ^  ^  ■  •  .  , ,  u 


resides  in  the 


waves  that  are  included  in  the  potential  <j>  but  are  ignored  in  the  zero-Froude- 

number  potential  .  For  the  present  problem  of  determining  the  wave  resistance, 

u  (1) 
it  naturally  seems  desirable  to  retain  the  waves  in  the  potential  f  As  a  matter 

0 

of  fact,  comparison  of  the  wave-resistance  curves  associated  with  the  approximation 
r^  and  the  low-Froude-number  approximation  rp3^  in  figure  A  show  that  the  waves  in 


the  potential 


have  a  remarkable  effect  upon  the  wave  resistance.  The  most  im¬ 


portant  difference  between  the  wave-resistance  curves  r^^  and  r»F^  is  a  large  phase 

(1)  ZF 

shift.  Furthermore,  the  curve  r  is  slightly  lower,  in  the  mean,  than  the  curve 

Vpv  ,  and  the  amplitude  of  the  oscillations  in  the  curve  r^  are  somewhat  less  than 

(1) 

those  in  the  curve  rJr  . 

1  (1) 

Figure  A  shows  that  the  first-order  slender-ship  wave-resistance  curve  r 

passes  through  the  five  brackets  of  experimental  values^  at  the  Froude  numbers 
F=  .18,  .20,  .22,  .24,  and  .266,  whereas  the  curves  associated  with  the  low-Froude- 
number  approximation  the  zeroth-order  approximation  r^°\  and  the  Michell 

approximation  r^  are  usually  outside  (well  outside  in  some  cases)  these  experimental 
brackets.  The  points  corresponding  to  the  average  of  numerical  results  for  the  GBMK 
low-speed  theory  clearly  are  in  phase  with  the  curve  r£p\  indeed  fairly  close  to  it 
(somewhat  below)  as  was  noted  previously,  and  these  points  also  are  usually  outside 
the  experimental  brackets.  The  results  corresponding  to  the  modification  of  the  GBMK 
theory  by  Nakatake,  Toshima  and  Yamazaki  (1979),  in  which  a  local  nonuniform  flow 
transformation  is  used  to  incorporate  the  effect  of  the  nonuniformity  of  the  double¬ 
hull  flow  upon  the  propagation  of  waves,  are  not  significantly  different  from  the 
GBMK  results,  and  in  fact  are  quite  close  to  the  curve  r^^  (except  for  values  of 
the  Froude  number  above  .24,  where  a  slight  phase  shift  may  be  observed). 

Also  shown  in  figure  A  is  the  wave-resistance  curve  corresponding  to  the 
simplified  low-Froude-number  asymptotic  form  of  the  GBMK  theory  proposed  by  Baba  (1979). 
This  resistance  curve  is  plausible  at  sufficiently  low  values  of  the  Froude  number, 
say  for  F<.22,  but  the  curve  is  well  above  the  experimental  brackets  at  the  Froude 


numbers  F=.24  and  F=.266.  A  notable  feature  of  the  Baba  resistance  curve  is  that 
its  oscillations  are  considerably  less  than  those  of  the  other  curves;  it  also 
appears  to  increase  more  rapidly  as  the  Froude  number  increases.  Figure  A  also 
shows  five  values  of  the  wave-resistance  coefficient  obtained  by  Kitazawa  and 
Kajitani  (1979)  on  the  basis  of  a  modified  form  of  the  GBMK  theory  in  which  the 
boundary  condition  at  the  hull  surface  is  satisfied  more  accurately.  Two  of  these 
values  are  fairly  close  to  the  first-order  slender-ship  curve  r^\  but  the  other 
three  values  are  much  below  it;  the  results  corresponding  to  the  two  highest  Froude 
numbers  are  well  below  the  experimental  results. 

Figure  A  finally  shows  the  numerical  results  obtained  by  Dawson  (1979).  These 
results  are  somewhat  below  the  first-order  slender-ship  wave-resistance  curve  r^\ 
but  clearly  appear  to  be  in  phase  with  it.  Dawson’s  result  at  the  Froude  number  .266 
is  somewhat  below  the  experimental  bracket,  and  the  result  at  F=.20  is  on  the  low 
side  of  the  experimental  bracket.  This  suggests  that  Dawson's  results  may  in  fact 
be  somewhat  low.  The  first-order  slender-ship  wave-resistance  curve  r^\  on  the 
other  hand,  i3  slightly  on  the  high  side  of  the  experimental  brackets,  perhaps  as 

G 

a  result  of  the  various  approximations  used  for  the  purpose  of  simplifying  the 
numerical  calculations.  The  approximation  r^  and  the  numerical  results  of  Dawson 
thus  appear  to  be  roughly  in  agreement. 

In  summary,  on  the  basis  of  the  results  shown  in  figure  A,  the  first-order 
slender-ship  wave-resistance  approximation  r^  appears  to  yield  results  significantly 
different  from  those  predicted  by  the  Gueve 1-Baba-Mar uo-Kayo  low-speed  theory  or  the 
several  above-mentioned  modifications  of  this  theory,  and  in  better  agreement  with 
experimental  results.  As  a  matter  of  fact,  among  the  five  resistance  curves  and 
the  four  sets  of  numerical  values  shown  in  figure  A,  only  the  first-order  slender- 
ship  wave-resistance  curve  r^  passes  through  all  the  five  experimental  brackets. 


a.  I  wish  to  thank  Mr.  Cheng-Yo  Chen  for  permitting  that  figure  A  be  included  in 
the  present  report.  This  figure  is  excerpted  from  Mr.  Chen's  Ph.D.  dissertation, 
currently  in  preparation. 

2  2 

b.  The  Wigley  hull  is  defined  by  the  equation  y=±.05(l-4x  ) (l-256z  ),  where 
-.05<x<.05  and  0>z>-.0625. 

c.  The  approximation  r^j^  shown  in  figure  A  does  not  include  the  nonlinear  free- 
surface  correction  integral  in  expression  (6.6)  for  the  Kochin  free-wave  amplitude 
function.  This  free-surface  integral  was  neglected  for  simplicity. 

d.  The  average  values  shown  in  figure  A  for  the  wave-resistance  coefficient  predicted 
by  the  Guevel-Baba-Maruo-Kayo  low-speed  theory  have  been  determined  from  the  numerical 
results  given  page  60  of  the  Proceedings  of  the  Workshop  on  Ship  Wave-Resistance 
Computations  (1979).  Specif icially,  the  following  values  were  used: 
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Froude  number: 

.16 

.18 

.20 

.22 

.24 

.266 

Kitazawa  &  Kajitani: 

_ 

_ 

.55 

.41 

1.02 

.69 

Maruo  &  Suzuki: 

.20 

.42 

.62 

.45 

1.00 

- 

Miyata  &  Kajitani: 

- 

- 

- 

- 

- 

.61 

Mori  (potential  flow) : 

.20 

.47 

.59 

.41 

1.03 

.62 

Nakatake  (Baba) : 

.22 

.48 

.57 

.40 

1.03 

.66 

Nakatake  (Guevel) : 

.24 

.45 

.56 

.40 

- 

.65 

Yim: 

- 

- 

- 

- 

1.04 

.65 

average: 

.22 

.46 

.58 

.41 

1.02 

.65 

e.  For  simplicity,  the  nonlinear  free-surface  term  K  in  expression  (9.3)  for 

the  first-order  slender-ship  approximation  to  the  Kocfiin  free-wave  amplitude  function 
was  ignored  in  the  approximation  r(^)  shown  in  figure  A.  The  near-field  potential 
in  expression  (9.1)  for  the  first-order  slender-ship  potential  <j>(l)  was  also 
neglected,  as  can  be  justified  for  low  Froude  numbers.  Finally,  calculations  were 
further  simplified  by  replacing  the  wave  potential  given  by  equation 

(9.1b)  by  its  Michell  thin-ship  approximation,  as  is  approximately  justified  by  the 
fairly  small  differences  that  may  be  observed  in  figure  A  between  the  zeroth-order 

slender-ship  wave-resistance  curve  r^®)  and  the  Michell  curve  rw. 

M 

f.  The  experimental  brackets  shown  in  figure  A  correspond  to  the  highest  and  lowest 
of  11  experimental  values  consisting  of  the  9  values  given  page  52  of  the  Proceedings 
of  the  Workshop  on  Ship  Wave-Resistance  Computations  (1979)  and  the  2  values  given  in 
Gadd  and  Hogben  (1965) .  The  one  negative  value  given  page  52  of  the  Proceedings  of 
the  Workshop  for  the  Froude  number  .18  was  ignored  however.  Furthermore,  the  11 
experimental  values  at  the  Froude  number  .266  have  been  corrected  for  effects  of 
sinkage  and  trim  by  using  both  multiplicative  and  additive  corrections  determined  from 
the  numerical  results  obtained  by  Gadd  (1979),  Guevel ,  Delhommeau,  and  Cordonnier 
(1979,  1980),  Dawson  (1979),  and  Daube  (1980)  for  the  Wigley  hull  held  fixed  (no 
sinkage  and  trim)  and  left  free  to  heave  and  trim. 


Comparison  of  experimental  and  numerical  results  for  the  Wigley  hull 
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2.  Formulation  of  the  problem 

The  problem  examined  In  this  study  is  that  of  predicting  the  wave  resis¬ 
tance  experienced  by  a  ship,  with  length  L  say,  in  steady  rectilinear  motion, 
with  speed  U,  at  the  free  surface  of  an  otherwise  quiescent  sea  of  large  depth 
and  horizontal  extent.  Water  is  supposed  to  be  homogeneous  and  incompressible, 
with  density  p.  Viscosity  effects  are  ignored,  and  irrotational  flow  is  assumed. 
Effects  of  surface  tension,  wavebreaking,  and  spray  formation  at  the  ship  bow 
are  also  neglected. 

The  flow  caused  by  the  ship  is  rendered  independent  of  time  by  observing 

it  from  a  moving  system  of  coordinates  attached  to  the  ship.  Specifically,  the 

Z  axis  is  chosen  vert  al  and  pointing  upwards,  with  the  mean  sea  surface  taken 

as  the  plane  Z=0.  Tue  X  axis  is  in  the  mean  sea  plane,  parallel  to  the  direction 

of  motion  of  the  ship,  and  pointing  towards  the  bow.  The  Y  axis  is  in  the  mean 

sea  plane,  orthogonal  to  the  X  and  Z  axes,  and  oriented  so  that  the  X,Y,Z  axes 

form  a  right-handed  Cartesian  system  of  coordinates. 

Flow  variables  are  made  nondimensional  with  respect  to  the  length  L  and  the 

speed  U  of  the  ship,  and  the  density  of  water  p.  The  nondimensional  coordinates 

— ► 

(x,y  ,z) :  (X,Y,Z) /L  and  velocity  potential  <j>  =  $/UL,  where  X  and  $  are  dimensional, 
thus  are  defined,  and  the  nondimensional  flow  velocity  is  given  by 

(3?/3x,3<j>/3y  ,3<i>/3z);  (3$/3X,3<5>/3Y,3<l’/3Z)/U. 

2 

The  velocity  potential  <{>  satisfies  the  Laplace  equation  V  4>=0  in  the  flow 
domain.  On  the  hull  surface,  the  potential  <j>  satisfies  the  condition  3$/3n=nx, 
where  is  the  x  component  of  the  unit  vector  n  normal  to  the  hull  surface  and 
pointing  into  the  water,  and  3<})/3n=V<}) •  n  is  the  derivative  of  <J>  in  the  direction 
of  the  vector  n.  On  the  sea  surface,  given  by  the  equation  z=F  ( 3 <J> /3x- 1  V$  |  / 2), 
the  potential  <}>  satisfies  the  condition 

3£/3z  +  F2(32'?/3x2-3|V-ti|2/3x+7<t>-v!v,f  |2/2)  =  0, 

1/9  9  9  9 

where  F  is  the  Froude  number  defined  as  F=U/(gL)  and  |V(f>|  h(3<J>/3x)  +  (3<j>/3y)  + 

2 

(3$/3z)  is  the  square  of  the  magnitude  of  the  fluid  velocity  vector  V<t>. 

The  position  of  the  sea  surface  is  evidently  not  known  a  priori.  The 
sea  surface  boundary  condition  will  then  be  transferred  to  the  mean  sea  surface 
/.=()  by  using  a  Taylor  series  expansion  about  the  plane  z=0,  in  the  usual  manner. 

We  may  then  obtain  the  mcan-KoM-.su r I  ace  boundary  condition 

H/3z  +  F2324>/3x?  -  F2q  ( )  +  (KfS*)  =  0, 


(2.1) 
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where  the  term  q(if)  is  defined  as 

q(<?)  =  3  j  V4,  |  2/3x  -  V4>-v[\7*|2/2  -  (3<fr /3x- 1  Vtf.  |  2/2) 3  (3 4. /3 z+F23 /3x2) /?z  ,  (2.2) 

4  3 

and  the  term  0(F  ^  )  represents  terms  that  are  at  least  of  fourth  power  in  the 
Froude  number  F  and  of  third  power  in  the  potential  t  or  its  derivatives,  and 
may  presumably  be  neglected. 

In  accordance  with  the  transfer  of  the  sea-surface  boundary  condition  to 
the  mean  sea  surface,  the  Laplace  equation  will  be  satisfied  in  the  mean  flow 
domain,  d  say,  bounded  upwards  by  the  plane  z=0,  or  more  precisely  by  the  portion, 

•  say,  of  the  plane  z=0  outside  the  ship  hull  surface;  we  thus  have 

"72:  =  0  in  the  mean  flow  domain  d.  (2.3) 

The  boundary  condition  on  the  hull  surface  will  similarly  be  applied  on  the  mean 
hull  surface,  h  say,  that  is  the  portion  of  the  ship  hull  located  below  the  plane 
z=0;  this  yields 

3; /in  =  n^  on  the  mean  hull  surface  h.  (2.4) 

A  major  difficulty  of  the  above  problem  stems  from  the  sea-surface  boundary 

2 

condition  (2.1,  2.2),  which  is  nonlinear.  However,  the  nonlinear  term  F  q($)  in 

2  2 

equation  (2.1)  is  of  order  F  tf>  ,  as  may  be  seen  from  equation  (2.2),  and  thus 

2  2  2 

can  be  presumed  to  be  small  in  comparison  with  the  linear  term  3<)>/3z+F  3  d>/3x 

2 

for  slender  ship  hull  forms.  The  nonlinear  term  F  q(<f>)  may  then  be  neglected 
in  a  first  approximation.  More  generally,  this  nonlinear  term  will  be  incorporated 
in  an  iterative  manner  by  expressing  the  sea-surface  boundary  condition  (2.1)  in 
the  form 

2  2  2  2 

3j/3z  +  F  3  y/3x  =  F  q(f)  on  the  mean  sea  surface  0,  (2.5) 

and  treating  the  term  F2q(f)  on  the  right  side  as  a  nonhomogeneous  term  for  the 

2  2  2 

linear  condition  .‘/•z+F  •  ‘  /  ■  x  =0. 

In  summary,  the  problem  of  flow  caused  by  steady  rectilinear  motion  of  a 
ship  in  a  calm  sea  is  formulated  in  this  study  as  the  potential-flow  problem 
defined  by  th-  I  '"lace  equation  (2.3)  subject  to  the  boundary  conditions  (2.4) 
and  (2.3)  on  th"  sh i p  hull  and  on  the  free  surface  of  the  sea,  respectively.  In 
addition,  the  usual  "radiation  condition",  specifying  that  waves  are  present 
only  behind  the  ship,  must  he  imposed  for  uniqueness  of  the  solution. 
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3.  The  Green  function 

The  above-defined  boundary-value  problem  will  be  solved  by  formulating 
an  integral  equation  for  the  velocity  potential  <j>  based  on  the  use  of  a  Green 
function  satisfying  the  linearized  sea-surface  boundary  condition.  Specifically, 
this  Green  function,  which  will  be  denoted  by  G(t,x;F2)  or  simply  by  G,  is  the 
linearized  velocity  potential  of  the  flow  created  at  point  £(£,n,£<0)  by  a  unit 
outflow  at  point  x(x,y,z<0),  stemming  from  a  submerged  source  if  z<0  or  from  a 
flux  across  the  mean  sea  surface  if  z=0. 

The  Green  function  satisfies  the  radiation  condit ion,  namely  that  there  are 
no  waves  for  £>x,  and  the  equations 


V  G  =  5  (x-£)6  (y-n)<$  (z-C)  in  z< 0, 

(3.1a) 

3G/3z  +  F2d2G/3x2  =0  on  z  =  0 , 

(3.1b) 

if  £<0,  and 

V2G  =0  in  z<0. 

(3.2a) 

3G/3z  +  F232G/3x2  =  -6  (x-O  6  (y-rj)  on  z  =  0, 

(3.2b) 

if  ?=0.  Equations  (3.1a,b)  are  well  known.  These  equations,  however,  are  valid 
only  if  C  is  strictly  negative,  and  equations  (3.2a,b)  are  proper  in  the  limiting 
case  C=0,  as  is  demonstrated  in  Noblesse  (1981). 

The  Green  function  may  be  expressed  in  the  form 

AirG(t,  x;F2)  =  -l/|t-x|  +  N(?,^;F2)  +W(£,x;F2),  (3.3) 

y  y  2  2  21/2 

where  the  first  term  on  the  right  side,  with  | £-x| = [ (£-x)  +(n-y)  +(C~z)  ]  , 

is  the  Green  function  for  potential  Flow  in  an  unbounded  fluid,  and  the  functions 
N(t,x;F2)  and  W(<t,x;F2)  stem  from  the  presence  of  the  sea  surface  and  represent 
a  nonoscillatory  near-field  disturbance  and  a  wavy  disturbance,  respectively. 

The  function  W,  representing  the  waves  behind  the  singularity  at  point  x,  is 
given  by  the  integral 

W  =  H(x-0(A/F2)£  Im  exp[(c+z)(l+t2)/F2+i{(e-x)+(n-y)t>(l+t2)Ji/F2]dt,  (3. A) 


where  H(x-0  is  the  Heaviside  unit-step  function,  which  takes  the  value  1  for 
x-?>0  and  0  for  x-£<0,  and  Im  represents  the  imaginary  part. 

The  function  N,  representing  a  nonoscillatory  near-field  disturbance,  is 
given  by  the  integral 


N  =  1  /  r '  + 


(2/nF2)  Tm 


exn(Z)E^(Z)dt  , 


(3.5) 
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where 

2  ?  21/2 
r'  =  [(C-x)z+(n-y)^+(r,+2)‘]1/z 

is  Che  distance  between  the  field  point  X  and  the  mirror  image  of  the  singularity 
point  x  with  respect  to  the  mean  sea  plane  z=0,  Z  is  the  complex  function  of  the 
real  variable  t  defined  by 

Z  5  [(c+z)(l-t2)1/2+(n-y)t+iU-x|  ](l-t2)1/2/F2, 

and  E^(Z)  is  the  exponential  integral,  which  is  defined  as  in  Abramowitz  and 
Stegun  (1965).  For  purposes  of  numerical  evaluation,  a  convenient  alternative 
expression  for  the  near-field  disturbance  N  is 

N  =  1/r '  -  (2/F2) [l-(c+z)/(r'+|C-x| )]  +  (2/ttF2)/^  Im[exp(Z)E1 (Z)+£n(Z)+y ]dt ,  (3.5a) 

where  y=0.577...is  Euler's  constant.  Expression  (3.5)  is  well  suited  for  obtain¬ 
ing  an  ascending  series  of  the  near-field  disturbance  N,  useful  for  evaluating 

'l 

N  for  small  values  of  r'/F  .  Indeed,  the  first  term  in  this  series  is  1/r',  and 
the  second  term  is  shown  in  expression  (3.5a).  The  third  term  in  the  series  is 
given  in  Noblesse  (1978),  where  the  exp-essions  for  VN  corresponding  to  express¬ 
ions  (3.5)  and  (3.5a)  may  also  be  found. 

An  interesting  alternative  form  of  expression  (3.5)  is 

N  =  -1/r'  +  (2/ttF2)[J  Im[exp(Z)E  (Z)-1/Z]dt.  (3.6) 

However,  this  expression  is  not  well  suited  for  purposes  of  numerical 

evaluation.  Neither  is  expression  (3.6)  suited  for  obtaining  an  asymptotic 

2 

expansion  of  the  near-field  disturbance  N  for  large  values  of  r'/F  .  A  comple¬ 
mentary  integral  representation  of  N  suited  for  obtaining  such  an  asymptotic 
expansion  is 

N  =  -1/r'  +  (2/tiF2)/“  Re[exp(Z')E  (Z')-1/Z']dt  +  (4/F2)J°°  Im  exp(Z')dt,  (3.6a) 

1  K~x|  / 1  n-y  I 

where  Re  represents  the  real  part,  and  Z'  is  the  complex  function  defined  by 
Z’  =  [  (r,+z)  (1+t2)  1/2+i  (  |  £-x|  -|  n-y  I  t)  )  (l+t^)1/2/F2. 

2 

Indeed,  an  asymptotic  expansion  for  large  values  of  r'/F  has  been  obtained  in 
Noblesse  (1975)  from  expression  (3.6a)  in  the  particular  case  when  n=y,  for 
which  the  last  integial  on  the  right  side  of  equation  (3.6a)  vanishes.  The 
first  term  in  this  asymptotic  expansion  is  (-1/r').  The  four  alternative 
integral  representations  of  the  near-field  disturbance  N  given  above  may  be 
found  in  Noblesse  (1981). 
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We  have 

N  -v  1/r '  as  r'/F2  ->  0  and  N  ^  -1/r'  as  r’  /F2  -»■  «>  .  (3.7a,b) 

Equation  (3.7b)  shows  that  far  behind  the  singularity  x,  as  (  +  -™  ,  the  near-field 
disturbance  -l/|£-x|+N  in  expression  (3.3)  for  the  Green  function  is  approximately 
equal  to  -2/(C  +nz+^^)-L/  .  This  is  negligible  in  comparison  with  the  wave  disturb¬ 
ance  W.  We  may  then  obtain 

G  (1/itF2)/^  Im  exp[(?+z)  (l+t2)/F2+i{(C-x)+(n-y)t)  (l+t2)1/2/F2]dt  as  r  (3.8 

This  result,  which  is  well  known  of  course,  will  be  used  further  on  in  this  study 
(in  section  6)  for  determining  the  amplitude  of  the  free  waves  far  behind  the  ship. 


\ 


»v. 


1  2 


4.  Basic  integral  identities  for  the  velocity  potential 

In  this  section,  basic  integral  identities  for  the  velocity  potential  are 
obtained  by  applying  a  classical  Green  identity  to  the  potential  <t>  =  $(x)  and  the 
Green  function  G2G(f,,x)  defined  above.  This  Green  identity  is 
o  2 

/,,  (<}>V~G-GV  4> )  dv  =  /  ,  (if  3G/3z-G34>/3z)dxdy  +  f,(G3<|>/3n-<J>3G/3n)  da 
a  *  o  J  n 

+  /,  (4>3G/3n-G8ij>/3n)  da,  (4.1) 

'  n 

CD 

where  d'  is  the  finite  domain  bounded  by  the  hull  surface  h,  the  mean  sea  plane 
z=0,  and  some  surface  h^  surrounding  the  hull  h,  as  is  shown  in  figure  1;  o'  is 
the  portion  of  the  plane  z=0  between  the  intersection  curves  c  and  c  of  the 

CD 

hull  surface  h  and  the  exterior  surface  h  ,  respectively,  with  the  plane  z=0.  On 
the  surfaces  h  aid  h  ,  we  have  3(j;/3n:iVd>*n  and  3G/3nEVG*n  where  n  is  the  unit  out- 

oo 

ward  normal  vector  to  h  or  h  ,  as  is  shown  in  figure  1.  Furthermore,  dv  and  da 

represent  the  differential  elements  of  volume  and  area  at  point  x  of  the  domain 

d'  and  the  surfaces  h  or  h  ,  respectively. 

By  expressing  the  integrand  <j>3G/3z-G3<(>/3z  of  the  first  integral  on  the  right 

2  2  2  2  2  2 

side  of  equation  (4.1)  in  the  form  ({’(SG/Sz+F  3  G/3x  )  -G(3<J>/3z+F  3  $/3x  )  + 

2 

F  3  (G3<}>/3x-<t>3G/3x) /3x,  and  by  using  the  relation 

j  ,3(G34>/3x-(J>3G/3x)/3xdxdy  =  /  (G3<j>/3x-<J>3G/3x)dy  +  /  (G3iJ>/3xH>3G/3x) dy , 

^OO 

we  may  obtain 

fd,<p\7ZGdv  -  f0,<P(  3G/3z+F232G/3x2)dxdy  =  /d,  GV2<pdv  -  /q  ,  G(3<t>/3z+F2324-/3x2)  dxdy 
+  /h(G3<j)/3n-({i3G/3n)da  +  F  /c  (G3<j>/3x-4>3G/3x)dy  +  1^ ,  (4.2) 

where  I  is  defined  by 

OO  J 

I  =  fu  (4>3G/3n-G3<f>/3n)da  +  F2/  (G3<J>/3x-<J>3G/3x)  dy . 

°°  '  n  J  c 

oo  oo 

It  can  be  shown  that  I  vanishes  as  the  exterior  surface  h  is  made  ever  larger. 
This  term  can  then  be  ignored  in  equation  (4.2),  and  the  finite  domain  d'  and 
region  o'  of  the  mean  sea  plane  may  be  replaced  by  the  unbounded  mean  flow  domain 
d  and  sea  surface  a,  respectively. 

By  expressing  the  potential  $  in  the  integrands  of  the  two  integrals  on 
the  left  side  of  equation  (4.2)  in  the  form  <f>  =<$>*+($-$*)  ,  where  <f>H<f>(x)  ns  was 
defined  previously,  and  4>A  represents  the  potential  at  point  t,  i.e.  <)>*-<l>(t),  we 
may  obtain 

JdtV2Gdv  -  /  4-( 3G/3z+F232G/3x2)dxdy  =  +  C*, 


(4.3) 
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where  C  and  C'  are  given  by 

C  =  /dV2G dv  -  /o(3G/3z+F232G/3x2)dxdy,  (4.4) 

C'  =  fd  (<P-<P  *)V2Gdv  -  /a  ($-<£*)  (  3G/3z+F23  2G/3x2)dxdy  . 

It  may  be  seen  from  equations  (3.1)  and  (3.2)  that  we  have  C';0  if  ,-*-0  as 
x-*£ ,  that  is  if  the  potential  is  continuous  everywhere  in  the  solution  domain 
d  and  on  its  boundary  o+h+c,  as  is  assumed  here.  Use  of  equation  (4.3),  with 
C'hO,  in  equation  (4.2)  then  yields 

=  J^GV  4>dv  -  ^  G ( 3<£/3z+F23 2<f> /3x  ) dxdy  +  (G3ij>/3n-<}>3G/3n)da 

+  F2/c(G3<fi/3x-4>3G/3x)dy,  (4.5) 

where  C  is  given  by  formula  (4.4). 

Use  of  equations  (3.1)  and  (3.2)  in  expression  (4.4)  for  C  shows  that  we 
have  C=1  if  the  point  £  is  in  the  mean  flow  domain  d  or  on  the  mean  sea  plane  a, 
but  strictly  outside  the  hull  surface  h,  whereas  we  have  C:0  if  £  is  strictly  in¬ 
side  h.  We  thus  have 

**  =  IdGv2t,dv  ~  !0  G(?<J>/3z+F2324>/3x2)dxdy  +  fh  (G3(J./3n-<j>3G/3n)da 

+  F  (  (G3<f>/3x-<j>3G/3x)  dy  (4.6) 

J  c 

for  X  in  d+a-h-c,  and 

2  2  2  2 

0  =  /  GV  «J>dv  -  /aG(34>/3z+F  3  <(>/3x  )dxdy  +  /h(G34>/3n-<j>3G/3n)da 

+  FZ[  (Gd<p/3x~<pdG/dx)dy  (4.6') 

for  £  in  d.+a.-h-c,  where  d^  and  represent  the  domain  and  the  portion  of  the 
plane  z=0,  respectively,  inside  the  hull  surface  h.  It  can  also  be  seen  from 
equations  (3.1)  and  (3.2)  that  we  have  C=l/2  if  the  point  X  is  exactly  on  the 
hull  surface  h  or  its  intersection  c  with  the  plane  z=0,  at  least  for  points  £ 
where  h+c  is  smooth;  more  generally,  the  value  of  4irC(or  2ttC)  at  a  point  t,  of  h 
(or  c)  is  equal  to  the  angle  at  which  d(or  a)  is  viewed  from  the  point  %.  We 
then  have 

<t>/2  =  f  ,GV2$dv  -  /  G(3<f>/3zl-F232(f>/3x2)dxdy  +  f  (G3<)>/3n-<t>9G/3n)da  (4.6") 

*  '  d  J  q  /  rj 

o 

+  F  /  (G3<t>/3x-tf>9G/9x)dy 
y  c 

for  l  exactly  on  smooth  h+c.  Equations  (A. 6)  are  well  known,  although  the  par- 

2  2  2  2 

ticular  case  of  these  equations  corresponding  to  V  (J)=0=3<^/^z+F  *  (f*/3x  is  usually 


given  in  the  literature,  and  these  equations  are  usually  obtained  from  the  Green 

identity  (4.1)  in  a  manner  different  from  that  shown  here. 

The  value  of  the  constant  C  on  the  left  side  of  equation  (4.5)  is  discontinuous 

across  the  hull  surface,  C  being  equal  to  1  outside  the  ship  and  to  0  inside,  as  is 

explicitly  indicated  in  equation  (4.6).  This  discontinuity  in  the  value  of  C 

evidently  is  accompanied  by  a  corresponding  discontinuity  on  the  right  side  of 

equation  (4.5).  Specifically,  the  latter  discontinuity  stems  from  the  integrals 

J  ^4*  3G/3nda  and  J^SG/Sxdy  representing  potentials  due  to  dipole  distributions  on 

the  hull  h  and  the  waterline  c,  respectively.  An  integral  identity  valid  for 
-► 

any  point  C  —  outside,  inside,  or  exactly  on  the  hull  surface  — —can  be  obtained 
by  eliminating  the  discontinuity  in  the  value  of  C  in  equation  (4.5).  This  can 
be  done  by  adding  the  term  C.<K  on  both  sides  of  equation  (4.5),  with  C.  given  bv 

C.  =  /d  V2C<JV  "  J  (^G/^z+K2a2G/'lx2)dxdy.  (4.7) 

i  i 

Equation  (4.5)  then  becomes 

=  /dG?%dv  -  JnG(3$/c>z+F232<}'/3x2)  dxdy  +  (G3<J>/3n-ij>3G/3n)  da 

+  F2/  (G3<f>/3x-<j>3G/3x)  dy  +  C.(j>.  (4.8) 

where  I  is  defined  as 

1  =  /d+d  v2cdv  ~  /0-k,  OG/3z+F232G/3x2)dxdy. 
i  i 

Equations  (3.1)  and  (3.2)  show  that  we  have  Isl  for  any  point  X  in  the  lower 
half  space  £<0.  By  using  the  divergence  theorem 

f,  V  Gdv  =  /,  ?G/3nda  +  /  3G/3zdxdy 
d .  'ii  'a, 

l  i 

and  the  relation 

/  32G/3x2dxdy  =  -/  3G/3xdy 
i 

in  equation  (4.7),  we  may  obtain 

C.  =  f.  3G/3nda  +  F 2f  3G/3xdy. 
l  '  n  '  c 

By  substituting  the  above  expression  for  into  equation  (4.8),  and  replacing 
I  by  1,  we  may  finally  obtain 

0*  =  fA GV2$dv  -  /oG(3i/3z+F232((!/3x2)dxdy  +  /  [G3<> /3n- (4>-<(. *)  3G/?n]da 

+  F2/c(G34>/3x-($-t*)3G/3x]dy  .  (4.9) 
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The  integral  identity  (4.9)  is  valid  for  any  point  £,  whether  outside,  inside, 
or  exactly  on  the  hull  surface  h.  This  new  identity  thus  is  essentially  equivalent 
to  the  set  of  the  three  classical  identities  (4.6),  (4.6'),  and  (4.6"),  which  are 
valid  exclusively  for  £  outside,  inside,  or  on  the  hull  surface  h,  respectively. 

As  a  matter  of  fact,  these  three  identities  can  be  derived  from  the  identity  (4.9) 
by  noting  that  we  have  /^3G/3nda  +  F  f^G/ 3xdy  =  0,  1,  or  1/2  for  £  outside,  inside, 
or  on  the  hull  surface  h,  respectively.  This  may  be  shown  by  writing  the  above 
expression  in  the  form 

/  3G/3nda  +  fQ  3G/3zdxdy  -  /  3C/3zdxdy  -  F 2f  32G/3x2dxdy  = 
i  i  l 

9  9  9  9 

/  V  Gdv  -  /  (3G/3Z+F  3  G/3x  )dxdy. 

i  i 

The  above-stated  result  and  the  three  integral  identities  (4.6),  (4.6'),  and  (4.6") 
then  readily  follow  from  equations  (3.1)  and  (3.2). 

An  interesting  feature  of  the  integral  identity  (4.9),  or  of  the  related  identities 
(4.6),  (4.6')  and  (4.6"),  is  the  appearance  of  a  waterline  integral  for  the  case  of 
a  sea-surface  piercing  hull  (no  waterline  integral  is  present  in  the  case  of  a  fully- 
submerged  body) .  A  simple  interpretation  of  the  waterline  integral  may  be  obtained 
by  considering  a  sea-surface  piercing  hull  as  the  "zero-submergence  limit"  of  the 
slightly-submerged  body  consisting  of  the  mean  hull  h  closed  by  a  horizontal  "lid"  t; 
specifically,  the  waterline  integral  can  then  be  shown  to  stem  from  the  effect  of  the 
lid  closing  the  slightly-submerged  body  h+£ .  Details  are  given  in  appendix  A.  The 
present  study  is  concerned  with  the  problem  of  potential  flow  about  a  ship,  that  is 
the  exterior  potential-flow  problem".  Integral  identities  corresponding  to  equations 
(4.6),  (4.6')  and  (4.6")  can  evidently  also  be  obtained  for  the  "interior  potential- 
flow  problem",  that  is  for  the  potential,  ^  say,  defined  in  the  interior  domain  d.. 
These  integral  identities,  as  well  as  identities  obtained  by  adding  the  integral 
identities  for  the  exterior  and  interior  problems,  are  listed  in  Appendix  B. 
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5 .  Equation  for  determining  the  velocity  potential 

Use  of  equations  (2.3),  (2.4),  and  (2.5)  in  the  fundamental  integral  identity 
(4.9)  yields 

=  /  Gn^da  -  /  (4'-<J *)  3G/3nda  +  F"-/  / Dx —  (<f  —4*^)  H;/''x]dy  -  F^/^Gq (>f> ) dxdy .  (5.1) 

Equation  (5.1)  is  the  basic  equation  obtained  in  this  study  for  determining  the 

velocity  potential  {■  on  the  mean  hull  surface  h+c  and  in  the  mean  fluid  domain 

d+c .  An  alternative  form  of  this  equation  will  now  be  given. 

Let  n  ,  n  ,  n  represent  the  x,  y,  z  components  of  the  previously-defined 
x  y  z  _>_ 

unit  outward  normal  vector  n  to  the  mean  hull  surface.  We  also  define  the  unit 
— ► 

vector  t  tangent  to  the  mean  waterline  c,  oriented  as  is  shown  in  figure  1;  the 
x,  v,  z  components  of  t  are  given  by  t  ,  t  ,  0.  If  the  unit  positive  vector  along 

->■  x  y 

the  x  axis  is  denoted  by  i,  we  then  have 

3f>/3x  =  V$.i  =  (n3^/?n+t34>/S£+nxt8t/^d)  .1  =  n^Siji /tn+t ^ / % l-n 3<t>/3d, 

where  ?f/;n  is  the  derivative  of  in  the  outward  normal  direction  n  to  h,  as  was 
defined  previously,  l  represents  the  arc  length  along  the  mean  waterline  c 
(oriented  as  is  shown  in  figure  1),  and  3<ji/9£  is  the  derivative  of  $  in  the  direc¬ 
tion  of  the  unit  tangent  vector  t  to  c;  finally,  3<|>/3d  represents  the  derivative 
of  $  in  the  direction  "d"  defined  by  the  unit  vector  nxt,  which  is  tangent  to  the 
hull  surface  and  pointing  downwards.  By  using  equation  (2.4),  we  may  then  obtain 

3 1 / 5 x  =  n2  +  tx  >;/•»£  -  nztv  ^7  'd. 

By  substituting  the  above  expression  for  3(J>/3x  into  the  waterline  integral  in 
equation  (5.1),  and  replacing  dy  by  t^dZ,  we  may  express  equation  (5.1)  in  the 
form 

4>(t)  =  1 5(t)  -  T ( t ; 4> )  •  (5.2) 

The  potential  ij)(t)  in  equation  (5.2)  is  defined  as 

Hi)  =  f.G n  da  +  F2/  Gn2t  dt,  (5.2a) 

n  x  c  x  y 

and  the  term  T(£;4>)  represents  the  transform  of  $  defined  below.  The  transform 
T(C;f>)  may  be  written  as  the  sum  of  a  linear  part,  L(t;<t>)  say,  and  the  nonlinear 
part  associated  with  the  nonlinear  free-surface  flux  q(<p).  Specifically,  we 
have 

T(£:f>)  =  L ( £ ; <f> )  +  F2/  Gq  0f> ) dxdy  , 

vT 


(5.2b) 


IS 


where  the  linear  transform  L(£,;<f)  is  given  by 


L  (£;<$)  =  (.(<fi-<t,JdG/dnda  +  F  2f  [  (<p-<p.)  dG/ dx-C  (t  3^/dl-n  t  3<p/dd)]t  dt .  (5.2c) 

'  n  *  *  q  *  x  z  y  v 

The  potential  $(0  In  equation  (5.2)  is  given  explicitly  in  terms  of  the  Froude 

number  and  the  geometry  (shape)  of  the  hull,  as  may  be  seen  from  equation  (5.2a), 

whereas  the  term  T(£;<}>)  is  of  course  not  known  a  priori. 

A  large  class  of  ships  operate  at  fairly  low  values  of  the  Froude  number.  It 
thus  is  interesting  to  briefly  examine  the  limiting  case  of  zero  Froude  number. 
Study  of  this  comparatively-simple  limiting  case  will  also  provide  useful  know¬ 
ledge  for  the  case  of  nonzero  Froude  number.  In  the  zero-Froude-number  limit, 
the  velocity  potential  41  becomes  the  zero-Froude-number  potential,^  say.  The 
potential  4>q  is  the  velocity  potential  of  the  flow  about  the  hull  surface  when 
the  sea  surface  is  replaced  by  the  rigid  plane  z=0,  as  may  be  seen  from  the 
zero-Froude-number  limit  of  the  sea-surface  condition  (2.5).  The  potential  4> 
thus  is  also  identical  to  the  potential  of  the  flow  of  an  unbounded  fluid  past 
the  hull  and  its  mirror  image  with  respect  to  the  plane  z=0  (the  double  hull) , 
and  indeed  is  often  referred  to  as  the  "double-hull  potential".  Specifically, 
the  zero-Froude-number  potential  satisfies  the  equations 


V  <|>0  =  0  in  d,  3<j)Q/3n  =  nx  on  h,  3<}>q/3z  =  0  on  a. 


(5.3a,b,c) 


corresponding  to  equations  (2,3),  (2.4),  and  (2.5),  respectively. 

The  equation  corresponding  to  equation  (5.2)  for  determining  the  zero-Froude- 
number  potential  ip q  takes  the  form 


t0<o  -  yo  -  yt-.y. 


(5.4) 


where  the  potential  <Pq(£),  corresponding  to  the  potential  $(£)  in  equation  (5.2), 


is  given  by 

=  /hGOnxda* 


'  1 


(5.4a) 


and  the  (linear)  transform  Lq(^;<}>q)  ,  corresponding  to  the  transform  T(C;4>),  is 
defined  as 


L0(?;*0)  =  /hU0(x)-*0(6]3G0/3nda.  (5.4b) 

In  equations  (5.4a,b),  Gq  is  the  zero-Froude-number  Green  function  given  by 
4tt G0(£,x)  =  -l/|f-x|  -  Hr', 
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where  r'  =  [ (£-x)  +  (n-y)  +  (C+z)  ]  as  was  defined  below  equation  (3.5). 

Equations  (5.4a)  and  (5.4b)  can  be  obtained  as  the  zero-Froude-number  limit 
of  equations  (5.2a)  and  (5.2b,c),  respectively.  Equations  (5.4)  may  also  be 
obtained  from  equations  (5.3a,b,c)  in  a  manner  analogous  to  that  used  for 
obtaining  equations  (5.2)  from  equations  (2.3),  (2.4),  and  (2.5):  specifically 
by  applying  the  Green  identity  (4.1)  to  the  zero-Froude-number  potential  <(>q 
and  Green  function  Gq ,  as  is  shown  in  detail  in  Noblesse  and  Triantafyllou 
(1980) . 

A  particular  result  of  the  latter  study  is  relevant  to  the  present  problem, 
namely  the  fact  that  we  have 

Lq (€  ; <Pq)  «  <PQ(t)  *  +Q(?),  (5.5) 

for  motion  of  a  slender  body  in  the  direction  of  its  major  dimension  (length), 
like  a  ship  in  forward  motion.  In  the  special  case  when  the  hull  surface  h  is 
the  lower  half  of  an  ellipsoid,  the  potential  $q(£)  actually  is  proportional 
to  the  exact  potential  ,  that  is  we  have  $q(£)  =  A$q(E;)  ,  and  the  value  of 

the  constant  of  proportionality  A  is  close  to  unity.  For  instance,  for  an 
ellipsoidal  hull  form  with  beam/length  and  draft/length  ratios  equal  to  .15 
and  .05,  respectively,  we  have  A - . 9 7 ;  in  other  words,  the  slender-ship  approxi¬ 
mation  $q(£)  is  smaller  than  the  exact  potential  <I>q(!)  by  about  3%.  It  is  also 
shown  in  Noblesse  and  Triantafyllou  (1980)  that  equation  (5.4)  may  be  solved 
efficiently  by  using  an  iterative  procedure  based  on  the  straightforward 
recurrence  relation 

/0k+1)(t)  -  V?)  -  VH<k>)'  '  '  <5.6) 

with  k>0  and  the  initial  (zeroth)  approximation  simply  taken  as  0, 

—  -V  ^  ( 1 )  ^ 

so  that  the  potential  i|L(C)  corresponds  to  the  first  approximation  in  the 

0  (k)  0 

sequence  of  slender-ship  approximations  <#>q  •  For  instance,  in  the  special 

case  of  the  above-defined  ellipsoidal  hull  form,  the  relative  error  associated 

(2)  -3 

with  the  second  iterative  approximation  is  approximately  10 

The  fact  that  we  have  iP(£)'''Vq(0  in  the  zero-Froude-number  limit  and 
$Q(t)r''4lQ(t)  in  the  slender-ship  limit  suggests  that  the  velocity  potential 
ijl(C)  defined  by  equation  (5.2a)  may  provide  a  useful  approximation  for 
common  slender  ships  operating  at  low  Froude  number.  Furthermore,  it  is 
interesting  to  note  that  the  potential  $(£)  becomes  identical  to  the 
classical  Michell  potential  in  the  thin-ship  limit.  More  precisely,  the 
Michell  potential,  <j>M(t)  say,  given  by 
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r\j(0  =  /hG(C,x,v=0,z)nxda, 

corresponds  to  the  first-order  approximation  in  a  th in-ship  expansion  of 
the  potential  1 0(C) ,  as  may  easily  be  seen. 

In  the  zero-Froude-number  and  thin-ship  limits  examined  above,  the  waterline 
integral  in  expression  (5.2a)  appears  to  be  negligible  in  comparison  with  the 
hull  integral.  However,  the  waterline  integral  may  be  expected  to  be  significant 
for  a  ship  form  with  a  blunt  bow  or /and  stern,  since  we  have  | n^  j  =1=  j  t^ |  at  a 
blunt  end.  The  waterline  integral  in  fact  is  of  utmost  importance,  as  has  been 
demonstrated  by  Eggers  (1980)  who  showed  that  the  velocity  field  associated  with 
the  hull  integral  in  expression  (5.2a)  is  not  continuous  along  the  mean  waterline 
c,  and  that  the  waterline  integral  is  precisely  needed  to  render  the  velocity 
continuous  along  c. 

The  potential  i l)(t)  defined  by  equation  (5.2a)  is  satisfactory  both  in  the 
thin-ship  limit  and  in  the  zero-Froude-number  slender-ship  limit,  and  the  velocity 
field  is  everywhere  continuous.  The  potential  $  therefore  seems  to  provide 
an  acceptable  slender-ship  approximation  to  the  potential  !.  More  generally, 
the  potential  ip  may  be  regarded  as  the  first-order  slender-ship  approximation 

in  the  sequence  of  slender-ship  approximations  <j/k^  defined  by  the  recurrence 
re  1 at  ion 

*(k+1)(?)  =  Ht)  -  T(t;0(k)),  (5.7) 


where  T(f;if>)  is  the  transform  defined  by  equations  (5.2b,c),  and  the  initial 
(zeroth)  approximation  is  simply  taken  as  <J>  ^  ?0. 

Equations  (5.2)  do  not  require  the  ship  to  have  port  and  starboard  symmetry; 
in  fact  the  hull  may  be  a  multiple  surface,  as  would  be  the  case  for  a  catamaran 
or  a  SWATH  ship.  In  the  most  common  case  of  a  single-hull  ship  with  port  and 
starboard  symmetry,  equations  (5.2)  may  be  simplified  somewhat  by  replacing 
the  mean  hull  h,  waterline  c,  and  sea  surface  o  by  their  corresponding  positive 
(port)  half  h+,  c+,  and  o+,  as  is  shown  in  figure  2,  and  the  Green  function 
G(E,x)  by  the  port-  and  starboard-symmetry  Green  function,  G(E,x)  say,  defined  as 

G(f , x)  =  G(C,x,y,z)  +  G ( f ,x , -y , z) . 

(1)  i 

Thus,  expression  (5.2a)  for  the  first-order  slender-ship  potential  <j>  - becomes 

1^^(\)  =  (,  Gn  da  +  F^f  Gn  t  <lf.  (5.8) 

x  x  v 

Numerical  evaluation  of  expressions  (.2a)  or  (5.8)  for  the  first-order  slender- 
ship  potential  will  hi-  examined  further  on  (in  section  9)  in  this  study. 
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6 .  The  Kochin  free-wave  amplitude  function  and 
the  Havelock  wave-resistance  formula 

The  wave  resistance  of  a  ship  can  be  determined  from  the  wave  pattern  far 
behind  the  ship,  as  is  well  known  since  Havelock  (1934).  The  velocity  potential 
<j>  in  the  "far  field"  (at  a  large  distance  away  from  the  ship)  may  be  determined 
in  terms  of  the  value  of  t  in  the  "near  field"  (on  the  mean  hull  surface  h+c  and 
on  the  mean  sea  surface  ;  in  the  vicinity  of  c)  by  means  of  the  equation 

4>(0  =  /,  (Gn  -<J>3G/3n)da  +  F  (  [G(n  +t  3<t>/3£-n  t  3<|>/3d)-<f>3G/3x]t  df  -  F“  f  Ga(0)dxdv.  (b.l) 
nx  '  c  xx  zy  v  c 

Equation  (6.1)  may  be  obtained  from  equations  (5.2)  by  replacing  the  terra 

'1/2 

i-i*  :  (x)-(f(t')  bv  $  in  equation  (5.2c)  since  we  have  <J>^->0  as  (£2+q2+t;2)  ,  and 

therefore  !  |  <<  |  <J>  |  for  f,  in  the  far  field  and  x  in  the  near  field.  In  other  words, 
equation  (6.1)  corresponds  to  the  far-field  limit  of  equation  (5.2).  However,  it 
is  interesting  that  equation  (6.1)  in  fact  is  valid  also  in  the  near  field.  More 
precisely,  equation  (6.1)  corresponds  to  equation  (4.6)  and  is  valid  everywhere 
in  the  domain  d+ T-h-c ,  that  is  for  f.  strictly  outside  the  hull  surface  h+c. 

If  we  are  only  interested  in  the  wave  pattern  far  behind  the  ship,  expression 
(6.1)  for  the  potential  MO  can  actually  be  greatly  simplified  by  replacing  the 
Green  function  G(\,x)  by  the  asymptotic  approximation  (3.8).  This  yields 

MO  ‘  (l/7r)Imj_vexp[  C(l+t?)/F2+i  (T+nt)  (l+t2)1^2/F2]K(t)dt  as  (6.2) 

where  the  function  K(t)  is  given  by 

_ n  n 

K(t)  =  F  ( En^— 3 E/ 3n) da  +  t^3cj>/ 3d)-<f>3E/3x]  t^d£  -  JoEq(<f)dxdv,  (6.3) 

with  the  function  ErE(x;t,F2)  defined  as 

E  =  exp[F  2(l+t2)^2{  (1+t2)  ^2z-i  (x+ty ) }]  .  (6.3a) 

The  equation  of  the  free  surface  far  behind  the  ship  is  given  by  C=F23<j>  (£  ,n ,  C=0)  /  3£ 
since  the  nonlinear  terms  in  the  Bernoulli  equation  may  be  neglected  at  a  suffi¬ 
ciently  large  distance  away  from  the  ship.  By  using  equation  (6.2)  we  may  then  obtain 

4(f»n)  ^  (l/Tr)Re/_OTexp[ i  (C+nt)  (1+t 2)''"/,2/F2 ]K(t )  (l+t2)^^2dt  as  £+-°°.  (6.4) 

Equations  (6.2)  and  (6.4)  express  the  potential  <ME)  and  the  equation  of  the 
free  surface  5(6, q)  far  behind  the  ship  in  terms  of  a  familiar  superposition  of 
elementary  plane  waves.  The  function  K(t)  essentially  gives  the  amplitude  of  the 
free  wave  component  at  angle  0=tan  ^ (  from  the  x  axis;  this  function  will  be  referred 
to  as  the  Kochin  free-wave  amplitude  function,  as  it  corresponds  to  a  particular 
case  of  the  function  introduced  by  Kochin  (1936)  for  determining  the  drag  and  lift 
acting  upon  a  ship  in  steady  rectilinear  motion.  The  nondimensional  wave  resistance. 
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r  say,  defined  as  r=R/oU2L2  where  R  is  the  dimensional  wave  drag,  may  be  directly 
obtained  from  the  Kochin  free-wave  amplitude  function  K(t)  by  means  of  the 
Havelock  (1934)  formula 

r  ?  R/pU2L2  =  (l/27T)J™jK(t)  I  2(l+t2)1/2dt.  (6.5) 

This  formula  may  be  obtained  by  performing  the  change  of  variable  t=tan6  in 
equation  (12)  page  118  of  Eggers,  Sharma,  and  Ward  (1967). 

Formula  (6.3)  for  the  Kochin  free-wave  amplitude  function  K(t)  may  be  ex¬ 
pressed  in  the  form 

K  (t )  =  F_2(l+t2)  Jhexp{F-2(l+t2)z}E[nx/(l+t2)-{nz-i(nx+tny)/(l+t2)1/2  .-:/F2]da 

+  f  E[n2+i(l+t2)1/24)/F2+t  3;J,/3d-n  t  3<j>/3d]t  dt  -  (  Eq(<j>)dxdv,  (6.6) 

'  c  x  x  zy  y  ' o 

2 

where  the  function  E-:E(x,y;t,F  )  is  given  by 

E  =  exp[-iF  2 (1+t 2) ^ ^2 (x+tv) ] .  (6.6a) 

Expression  (6.6)  for  the  Kochin  function  K(t)  can  be  evaluated  numerically  without 
difficulty  in  principle,  given  the  potential  <f>  on  the  mean  hull,  the  tangential 
velocity  components  34>/3d  and  3<(>/3d  at  the  mean  waterline  c,  and  the  nonlinear 
free-surface  flux  q  ( 4>)  on  a.  In  practice,  however,  the  rapid  oscillations  of  the 
exponential  term  E  for  small  values  of  the  Froude  number  F  or/and  large  values  of  t 
may  cause  a  loss  of  accuracy.  Furthermore,  some  terms  in  the  integrands  of  the  hull 
and  waterline  integrals  cancel  out  one  another,  as  is  shown  below.  This  requires 
that  comparable  methods  be.  used  for  numerically  evaluating  the  hull  and  waterline 
integrals.  For  instance,  it  may  be  useful  to  partially  group  the  two  integrals  in 
the  manner  shown  below,  or  in  some  similar  fashion. 

Let  the  mean  waterline  c  be  defined  by  the  parametric  equations 

x  =  afd)  and  y  =  b(d),  (6.7a,b) 

where  d  is  the  arc  length  along  c  as  was  defined  previously.  We  then  have 
t  =a'=da/d d  and  t  =b  '^db/dt.  We  define  the  "tangent  hull,"  hC  say,  tangent  to 

x  y  t 

the  mean  hull  h  along  the  mean  waterline  c.  Specifically,  the  tangent  hull  h  is 
defined  by  the  equations 

x  =  a (d)  +  za (d)  ,  y  =  b(d)  +  zB(d)  ,  0  •>  z  >_  -6(d)  ,  (6.8a,b,c) 

where  a(d)  and  B(d)  are  the  partial  derivatives  3x/3z  and  3y/3z,  respectively,  at 

the  waterline  c.  On  the  above-defined  tangent  hull  h  ,  we  have  nda=(3x/3zx3x/9d)d('dz. 

We  may  then  obtain  n  da=-(  b '+z8 ')d('dz,  n  da=(a'+zu')dldz  ,  and  n  da={ab  '-8a  '+z  (aB  ~-ba')  } 

x  y  -v  z  t 

dddz.  At  the  waterline  c,  the  unit  outward  normal  vector  n  to  h  is  given  by 
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2  1/2 

nx=-cb',  n^ca',  and  nz=c  (ab '-8a ')  where  c=l/ [l+(ab '-Ba ')  ]  .  If  the  potential 

:•  at  the  point  x=a+za,  y=b+zp,  z  is  taken  equal  to  the  potential  at  the  corres¬ 
ponding  point  x=a,  v=b,  z=0  on  c,  the  hull  integral  in  expression  (6.6)  can  be 
expressed  as  the  sum  of  an  integral  over  the  tangent  hull  ht  plus  an  integral  over 
the  hull  h+ht .  The  surface  integral  over  the  tangent  hull  h1”  can  be  partially 
integrated  analytically  in  the  form  of  a  line  integral  along  the  mean  waterline  c, 
which  can  be  combined  with  the  waterline  integral  in  expression  (6.6).  We  mav 
then  express  the  Kochin  function  K(t)  in  the  form 

K  (t )  =  /cE[c2(b')3+i(l+t2)1/2b'<|>/F2+a'b'9<}>/9.e+c(b')2(Ba'--ab')9<J>/9d 
-y(l-e)  (b  "/  (l+t2)+i  (b  '-ta')  <f>/F2(l+t 2)  1/2+(ctb  "-6a  ')  <*/F2} 

-y{6e-F  y(l-e)/(l+t  )}{B'/(l+t  )+i(6'-ta'H/F  (1+t  )  ^+(oi6'-6u')<>/F  }]d£ 

+  F_2  ( 1+t 2)  j  exp{F~2  ( 1+t 2)  z  }E[n  /(l+t2)-{n  -i(n  +tn  )  /  (l+t2)1/2  !(j>/F2]da 
h+hC  X  z  x  v 

-  !  Eq  (i)dxdy,  (6.9) 

Cf 

where  we  have 

E  =  exp [-iF_2 (1+t 2) 1/2 (x+ty) ] , 

as  was  defined  previously  in  equation  (6.6a),  and 

c  =  1/ [l+(ob  "-Ba') 2 ] 1/2  ,  y  =  (l+i(a+t6)/(l+t2)1/2]/[l+(a+tB)2/(l+t2)]  ,  and 
e  =  exp[-F~2(l+t2)6{l-i(a+t6)/(l+t2)1/2}). 

In  the  integral  over  the  tangent  hull  h1"  in  expression  (6.9),  the  potential  <J>  must 
be  taken  equal  to  the  potential  at  the  mean  waterline  c,  as  was  explained  previously; 
furthermore,  the  unit  normal  vector  n  must  be  taken  to  be  pointing  inwards,  rather 
than  outwards. 

The  integrand  of  the  integral  over  the  hull  h+ht  in  expression  (6.9)  vanishes 

as  z-*0  since  the  hull  h  and  the  associated  tangent  hull  h*1  are  tangent  along  the 

mean  waterline  c  and  have  opposite  unit  normal  vectors  n,  as  was  noted  above. 

-2  2 

Furthermore,  the  exponential  factor  expiF  (1+t  )z)  renders  the  integrand  small 

2  2 

for  negative  values  of  z,  especially  if  we  have  (1+t  )/F  >>1.  The  modified  hull 

t  2  2 

integral  over  h+h  in  expression  (6.9)  therefore  vanishes  as  (1+t  )d/F  ->»,  where 

d  is  of  the  order  of  the  draft  of  the  ship,  and  this  surface  integral  can  be 

neglected  in  comparison  with  the  modified  waterline  integral  in  expression  (6.9) 

2  2 

for  sufficiently  large  values  of  the  parameter  (1+t  )d/F  . 

Typical  values  of  the  Froude  number  and  of  the  draft  for  slow-speed  ships  may 

2  2 

be  taken  as  F=.15  and  d=.05.  This  yields  d/F‘"=;2.22  and  exp(-d/F  )  =  .ll,  which 
indicates  that  the  hull  integral  in  expression  (6.9)  may  be  neglected  in  comparison 
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with  the  waterline  integral,  as  a  first  approximation.  Furthermore,  the  slopes 
3x/3z=a  and  3y/3z=B  of  the  mean  hull  h  at  the  waterline  c  are  small  for  common 
low-speed  ships,  so  that  the  last  terms  (those  involving  a'  and  6")  in  the 
waterline  integral  in  expression  (6.9)  may  be  neglected.  Finally,  the  term 
c(b')  (6a'-ab')34>/3d  also  appears  to  be  negligible  for  common  low-speed  ships, 
partly  due  to  the  fact  that  8  and  a  are  usually  small  as  was  just  noted,  and  that 
we  have  3<j>/3d=~34>/3z=0(F  ).  Under  the  above  approximations,  we  may  then  obtain 
the  following  approximate  expression  for  the  Kochin  function  K(t): 

K(t)  -  |cE[{c2(b')2-Y(l-e)/(l+t2)+a'^/3f}b' 

+Y(l_e){i(l+t2)1/2bVY(l-e)-i(b'-taO/(l+t2)1/2+Ba'-ab'H/F2]d>e,  (6.10) 

where  the  nonlinear  free-surface  correction  integral  has  been  ignored.  A  main 
recommendation  of  the  above  approximate  expression  for  the  Kochin  function  is 
its  relative  simplicity.  Indeed,  equation  (6.10)  expresses  the  Kochin  function 
as  a  single  integral  on  the  mean  waterline  c,  and  only  requires  the  value  of  the 
potential  <f>  on  c. 

2  2 

The  smaller  the  values  of  F  /d(l+t  ),  |a|  and  |b|,  the  more  accurate  the 
approximate  expression  (6.10)  for  the  Kochin  function  K(t)  will  be.  In  particular, 
the  tangent  hull  h11  is  identical  to  the  hull  h  in  the  case  of  a  hull  in  the  form 
of  a  cylinder  orthogonal  to  the  plane  z=0,  and  the  hull  integral  in  expression 
(6.9)  vanishes  exactly  if  the  potential  4>(x)  is  taken  as  the  zero-Froude-number 
potential  ^(x.y).  By  putting  ot=0,  6=0,  c=l,  y=l ,  and  e=0  in  equation  (6.10), 
we  may  then  obtain 

K(t)  =  |cE[{(b02-l/(l+t2)+a'34)/3£}b"+it(a>+tb')4>/f‘2(l+t2)1/2]d£.  (6.11) 

21/2  2  2  21/2 

The  terms  i(l+t  )  t  $/F  and  in  $/F  (1+t  )  in  the  waterline  and  the  hull 

y  x 

integrals,  respectively,  in  expression  (6.6)  have  been  combined  into  the  term 

2  22  1/2 

it  b'<j>/F  (1+t  )  in  expression  (6.11).  The  fact  that  the  latter  term  vanishes 

as  t-*0  evidently  indicates  that  the  abovementioned  terms  in  the  waterline  and 

the  hull  integrals  in  expression  (6.6)  cancel  out  one  another  as  t-*D.  For  large 

2  1/2  2 

values  of  (1+t  )  /F  the  major  contribution  to  the  integral  (6.11)  stems  from 


points  of  stationary  phase,  if  any.  Points  of  stationary  phase  are  defined  by  the 

2  2  1/2 

relation  dx/d£+tdy/d£  :  a'+tb  '-0.  The  term  i  t  (a'+tb '  ).|>/F  (1+t  )  in  equation  (6.11), 

21/2  2  2  21/2 

which  stems  from  the  terms  i(l+t  )  t  <Ji/F  and  i  (n  +tn  )<j>/F  (1+t  )  in  the 

y  x  y 

waterline  and  the  hull  integrals  in  equation  (6.6),  thus  vanishes  at  a  point  of 
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3  2 

stationary  phase.  The  terms  00  and  b'/(l+t  )  in  expression  (6.11),  corres- 

2  2 

ponding  to  the  terms  n  t  and  n  / (1+t  )  in  the  waterline  and  the  hull  integrals 
'  x  y  x 

in  expression  (6.6),  may  also  be  shown  to  cancel  out  at  the  point  of  stationary 
phase  t=-a' /b' .  These  cancellations  of  terms  in  the  waterline  integral  and  in 
the  hull  integral  in  expression  (6.6),  for  t=0  and  at  points  of  stationary  phase, 
indicate  that  partial  combination  of  the  hull  and  waterline  integrals,  in  the 
manner  shown  above  or  in  some  similar  fashion,  may  be  useful.  The  cancellations 
also  demonstrate  the  importance  of  the  waterline  integral. 

In  the  common  case  of  a  single-hull  ship  with  port  and  starboard  svrametrv, 
the  Kochin  free-wave  amplitude  function  K(t)  is  an  even  function  of  t,  and  the 


Havelock  formula  (6.5)  becomes 
r  =  (1/ir)  J"|K(t)  |2(l+t2)1/2dt. 


(6.12) 


Furthermore,  expressions  (6.6)  and  (6.9)  for  the  Kochin  function  may  be  expressed 
in  terms  of  integrals  over  the  positive  (port)  halves  c+,  h+,  o+  of  the  mean 
waterline  c,  hull  h,  and  sea  surface  o. 

The  Havelock  wave-resistance  formulas  (6.5)  or  (6.12),  expressions  (6.6)  or 
(6.9)  for  the  Kochin  free-wave  amplitude  function  K(t),  and  equation  (5.2)  for 
determining  the  velocity  potential  the  "near  field"  of  the  ship  form  a  complete 

set  of  equations  for  evaluating  the  wave  resistance  of  a  ship  in  steady  rectilinear 
motion  in  a  calm  sea,  given  the  speed  of  the  ship  and  the  shape  and  position  of  its 
hull.  The  position  of  the  hull  is  not  known  precisely  beforehand,  due  to  the  sinkage 
and  trim  experienced  by  the  ship.  Equations  for  determining  the  sinkage  and  trim, 
which  requires  evaluation  of  the  hydrodynamic  lift  and  moment  exerted  upon  the 
ship,  must  then  be  added  to  the  above  set  of  equations.  However,  equations  for 
determining  the  hydrodynamic  lift  and  moment  and  the  resulting  sinkage  and  trim  will 
not  be  considered  explicitly  in  this  study,  which  is  concerned  with  the  basic  hydro- 
dynamical  problem  of  predicting  the  flow  caused  by  a  ship  and  the  ship  wave  resistance, 
given  its  speed  and  the  shape  and  position  of  its  hull. 

Equation  (5.2)  for  determining  the  velocity  potential  <j>  may  be  solved  by  using 

an  iterative  procedure,  as  was  discussed  in  the  previous  section.  Specifically,  a 

sequence  of  slender-ship  iterative  approximations  <j>^=  0,  i?,  may  be 

defined,  as  is  indicated  by  the  recurrence  relation  (5.7).  Corresponding  sequences 

of  slender-ship  approximations  K  K^\...  and  r^°\  r^\...  to  the  Kochin  free- 

wave  amplitude  function  K  and  the  wave  resistance  r  can  readily  be  defined,  by  simply 

(k)  (k) 

replacing  <p  by  <J>  and  K  by  K  ,  with  k  =  0,1,...,  in  equations  (6.6)  and  (6.5). 

The  zeroth-,  first-,  and  second-order  slender-ship  wave-resistance  approximations  are 
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given  and  examined  in  some  detail  in  the  following  sections.  The  classical 
Michell  (1898)  and  Hogner  (1932)  wave-resistance  formulas  and  the  low-Froude- 
number  wave-resistance  formulas  proposed  by  Guevel ,  Vaussy ,  and  Kobus  (1974), 
Baba  (1976),  Maruo  (1977),  and  Kayo  (1978),  which  can  be  related  to  the  zeroth- 
and  first-order  approximations  r^  and  r^  \  respectively,  will  also  be  con¬ 
sidered. 
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7 .  The  zeroth-order  slender-ship  wave-resistance  approximation 

In  this  section,  the  zeroth-order  slender-ship  wave-resistance  approxima¬ 
tion  r  0  corresponding  to  the  zeroth  approximation  <f>^^=0  is  considered.  By 
putting  <j>H0  into  expression  (6.6),  we  may  obtain  the  following  expression  for 
the  zeroth-order  slender-ship  approximation  to  the  Kochin  free-wave  amplitude 
function: 

K(0)(t)=/  En2t  d t  +  F-2L  exp{F-2(l+t2)z}En  da,  (7.1) 

•'cxy  n  x 

2 

where  E  =  E  (x,  y ;  t  ,F“)  is  the  exponential  function  defined  by  formula  (6.6a). 

An  alternative  expression  for  K^^(t)  convenient  for  computational  purposes 
is  given  by  expression  (6.9): 

K(0)(t)  =  /  E[c2(b')2-Y(l-e)/(l+t2)]b"d£  +  F-2/  exp{ F~2 (1+t 2) z }En  da,  (7.1a) 
C  h+h  X 

where  h  is  the  tangent  hull  associated  to  h  and  defined  by  equations  (6.8), 
and  the  functions  c,  y,  and  e  are  defined  below  equation  (6.9).  The  waterline 
integral  in  formula  (7.1a)  is  expressed  in  terms  of  the  shape  of  the  mean  water¬ 
line  c,  defined  by  the  functions  a(£)=x  and  b(£)=y,  the  slopes  a(l)=dxldz 
and  8(£)H3y/3z  of  the  mean  hull  h  along  the  waterline  c,  and  the  "equivalent-draft 
function"  <$(£)  related  to  the  local  draft  and  shape  of  the  framelines  of  the 
ship.  The  hull  integral,  on  the  other  hand,  depends  on  the  precise  shape  of 
the  "lower  hull",  at  some  depth  below  the  mean  waterline.  This  integral  is 
negligible  in  comparison  with  the  waterline  integral  for  sufficiently  large 
values  of  (l+t2)d/F2. 

The  above-defined  zeroth  approximation,  corresponding  to  the  trivial 

slender-ship  approximation  f’iO,  may  seem  extremely  crude.  It  will  however 

be  noted  that  the  approximation  (7.1)  corresponds  to  the  free  waves  associated 

with  the  first-order  slender-ship  potential  given  by  equation  (5.2a). 

Further  support  for  the  zeroth  approximation  is  provided  by  considering 

expression  (6.11)  for  the  Kochin  function  K(t),  obtained  in  the  particular 

2 

case  when  a-MD,  B"K),  and  F  /d-K),  for  the  limiting  case  when  t-K),  corresponding 

to  the  transverse  waves  in  the  spectrum  of  free  waves  following  the  ship.  By 

2  2 

using  the  relation  (b')  -l-(a')  in  expression  (6.11),  we  may  obtain 
K(t=0)  -Jcexp(-ix/F2)  (a  >-3iJi/3£)a'b  'dt  as  a-+0,  8->0,  and  F2/d-K).  (7.2a) 

In  the  same  limiting  case,  expression  (7.1a)  becomes 
K^(t=0)  <v  -Jcexp(-ix/F2)  (a')2b"d£  as  a-K),  0->O,  and  F2/d-K). 


(7.2b) 
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For  a  slender  ship,  with  beam/length  and  draft/length  ratios  of  order  e  sav, 

2 

we  have  d$/Zl=0(c  )«a'=0(l).  At  the  bow  and  stern  of  a  blunt-ended  ship 
form,  both  a'  and  3<j>/3-£  vanish;  however,  the  result  d<b/d&«a'  remains  valid, 
as  may  be  shown  by  considering  the  particular  case  of  flow  about  an  ellipsoid, 
for  instance.  It  may  then  be  seen  from  equations  (7.2a,b)  that  for  common 
low-speed  ships,  for  which  a,  8,  and  F2/d  are  small,  we  have  K^^(t)-K(t)  for 
small  values  of  t.  Finally,  it  is  interesting  that  the  classical  wave-resistance 
formulas  of  Hogner  (1932)  and  Michell  (1898)  correspond  to  particular  limiting 
cases  of  the  zeroth-order  slender-ship  approximation  K  ^  defined  by  formula  (7.1), 
as  is  discussed  in  some  detail  below. 

For  a  ship  with  a  fine  waterline,  that  is,  if  the  angle  between  the  mean 

2 

waterline  c  and  the  x  axis  is  sufficiently  small,  we  have  nx I t^ | << | | <<1 ,  and 
the  waterline  integral  in  formula  (7.1)  may  be  neglected  in  comparison  with  the 
hull  integral.  The  latter  integral  in  fact  is  identical  to  the  approximation, 

Kjj  say,  proposed  by  Hogner  (1932),  that  is  we  have 


^(t)  =  F  2/hexp{F  Z(l+tZ)z}exp{-iF  Z(l+tZ)X/ Z(x+ty) In^da. 


-,-2 


-2, 


.2,1/2, 


(7.3) 


The  Hogner  approximation  may  thus  be  obtained  as  the  fine-waterline  limit  of 
the  zeroth  approximation  K  ^  . 

In  the  case  of  a  single-hull  ship  with  mean  hull  h  defined  by  the  equation 
y=±b(x,z),  the  surface  integral- in  equation  (7.3)  can  be  transformed  into  the 
double  integral 


K^t)  =  2//h  exp  [F  2(l+t2)'L^2{  (l+t2)1/,2z-ix}]cos[F  2(l+t2) ^2tb  (x,z)  ] 3b (x , z)  /  -xdxdz  , 

y 

where  h  is  the  projection  of  h  onto  the  ship  centerplane  y=0.  If  the  ship  is 

y  -221/2, 

thin,  that  is  if  b(x,z)  is  sufficiently  small  that  the  term  cos[F  (1+t  )  “tb(x,z)J 
may  be  approximated  bv  1,  tiie  Hogner  approximation  .becomes  the  famous  Michell 
thin-ship  approximation,  say,  first  obtained  by  Michell  in  1898  and  given  bv 

KM(t)  =  2//hexp[F~2(l+t")1/,“{  (l+t2)1/,2z-ix}]3b(x,z)/3xdxdz.  (7.4) 

y 


More  precisely,  the  Michell  approximation  (7.4)  corresponds  to  the  first-order 
approximation  in  a  consistent  thin-ship  expansion  of  the  zeroth  approximation  K 
If  the  equation  of  the  mean  hull  h  is  expressed  in  the  form  z=-d(x,y),  the 
surface  integral  in  equation  (7.3)  can  be  transformed  into  the  double  integral 

K.„(t)  -  ff.  exp { -F  (l+t“)d(x,v) }exp{-iF  “(1+t  )  (x+t v ) > 3d (x ,y ) / Sxdxdv , 

H  J  n 

z 


(0) 
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where  h  is  the  projection  of  h  onto  the  mean  sea  plane  z=0.  If  the  ship  is 

Z  -2  2 
fiat,  that  is  if  d(x,v)  is  sufficiently  small  that  the  term  exp{-F  (1+t  )d(x,y) > 

mav  be  approximated  bv  1,  we  may  obtain  the  flat-ship  approximation,  say, 

ri 

proposed  by  Hogner  (1932)  and  given  by 

K^(t)  =  //^expi-iF  2(l+t2)1//2(x+ty)  }Sd(x,y)/9xdxdy.  (7.5) 

z 

The  Hogner  approximation  (7.3)  in  fact  was  obtained  by  Hogner  (1932)  as  an 
interpolation  formula  between  the  flat-ship  approximation  (7.5)  and  the  Michell 
thin-ship  approximation  (7.4). 

-2  2  1/2 

If  the  ship  is  both  thin  and  flat,  that  is  if  the  term  exp{-iF  (1+t  )  ty} 

-2  2 

in  formula  (7.5)  for  the  Hogner  flat-ship  approximation  and  the  term  exp{F  (1+t  )z) 
in  formula  (7.4)  for  the  Michell  thin-ship  approximation  may  be  approximated  by 
1,  both  of  these  approximations  yield  the  Maruo-Tuck-Vossers  slender-ship  approxi¬ 
mation,  Kmtv,  given  by 

SlTV  =  J3°ernexp{-iF"2(1+t2)1/2x}A^x)dX>  (7*6) 

where  A' (x) ^dA(x) /dx  and  A(x)  is  the  cross-sectional  area  of  the  mean  hull  h. 

The  slender-ship  approximation  (7.6)  was  obtained  by  Maruo  (1962),  Tuck  (1964), 

and  Vossers  (1962)  by  using  the  method  of  matched  asymptotic  expansions. 

-2  2  1/2 

The  thin-ship  approximation  exp{-iF  (1+t  )  tyj'vl  and  the  flat-ship 
-2  2 

approximation  exp{F  (1+t  )zhl  are  obviously  not  uniformly  valid  with  respect 
to  t  and  F.  As  a  matter  of  fact,  differences  between  the  approximations 
Kh,  and  can  be  quite  large.  This  may  be  seen  clearly  in  the  case 

of  a  simple  planar  hull  for  which  the  foregoing  approximations  can  be  evaluated 
analytically.  We  consider  the  hull  form  defined  by  the  equation 

y=t (b/2) (1-2 | x | +z/d) ,  where  0>z>-d ( 1-2 | x | ) ,  -l/2<x<l/2, 


and  b  and  d  are  constants  representing  the  beam/length  ratio  and  the  draft/length 

2  2 

ratio,  respectively.  The  nondimens ional  wave  resistance  r-R/pU  L  may  be  expressed 
in  the  form 


r  =  (16b2d2/Tr)f°k2(t)  ( 1+t 2T1  /2dt, 

2  1  /2 

whore  the  function  k(t)  i(l+t  )  K(t)/4bd  is  given  below, 
no tat i on 

(= ( 1+t 2 ) 1 7  2/2F2 ,  M>t(l+l2)1/2/2F2,  5=d(l+t2)/F2, 


For  shortness,  the 


will  be  used. 


so 


The  function,  k^  say,  corresponding  to  the  waterline  integral  in  the 
zeroth-order  approximation  (7.1)  is  given  by 


k...  =  4b-dF"  (cos.i-coss)  /  (  l-b“t“)  (b2+4d2+4b2d2)  . 


The  function,  k^ ,  corresponding  to  the  hull  integral  in  expression  (7.1)  and 
the  Hogner  approximation  (7.3)  is  given  by 


[b  t‘"+4d~(l+L  )  ]k  =  b“t  |  (s  inh)  /  fi-  (s  inn)  /a]  /  (l-b2t2)  -  4dF~  ( 1-cos  <)/[  l+4d~  ( 1  +  t  "I 


+  4d~  (1+t  “)  [  ( 1- 1 /«.  yp  '  I  /  ' -(sin  x) /a] /[  l+4d2  (1+t2)  ]  +  4dF“(c  osB-cosa)  /  1 1  -h“t  1  . 


(0) 


>rresponding  to  the  zeroth-order  approximation  (7.1)  is 


The  function,  k 

given  by  k^  ^=k^+k^.  By  adding  the  above  expressions  for  kw  and  k^ ,  we  may 
then  obtain 


">  >  1  i  (01  2  7  7  7  7  77 

[b*-t^+4d“(l+t‘)]kV  ;  =  b  t  [(sinp)/e-(sina)/a]/(l-b  t  )  -  4dF^  (1-cosct) / [l+4d  (l+t“) ! 


+  4d“ ( 1+t 2 ) [ (1-1/ex p: ) / 3- (sina) /a] / [ 1+4 d2 (1+t2) ]  +  4dF2(b2+4d2) (cosF-cos  •) / (b2+4d2+4b2d^ 


Th 


e  function,  k^,  corresponding  to  the  Michell  approximation  (7.4)  is  given  by 


K 


( 1-1 /expo)  /i-(sina)  /  t+4dF2  (1-cosa)  ]  /  [l+4d2  (1+t ")  ] , 


(0) 


which  is  identical  to  the  limit  b=0  of  the  functions  k^  and  kvu/.  The  function, 
k(j,  corresponding  to  the  Hogner  flat-ship  approximation  (7.5)  is  given  by 

k^  =  [ (sinB) /B-(sina) /u] / (l-b2t2) , 

which  is  identical  to  the  limit  d=0  of  the  functions  k^  and  k  ^ .  Finally,  the 
function,  k^.^,  corresponding  to  the  Maruo-Tuck-Vossers  approximation  (7.6)  is 
given  hy 


kMTV  =  1- (s ina) /a , 


which  is  identical  to  the  limits  b=0  and  d=0  of  the  functions  and  k^,  respect ivelv . 
Differences  between  the  above  functions  are  most  striking  and  important  for 


large  values  of  t.  In  the  limit  t the  waterline  integral  k  and  the  hull  integral 

w 


(Hogner  approximation)  k  become 

H 


k^  *  F2  (cosB-cosa)  /  (d+b^4d+b  2d)  t 2  as  t  a-  <*>, 


k^  v  F2  [  (cosa) /d-2  (bsinB+2dcosB)  /  (b2+4d“')  ] /b2t^  as  t 


We 


thus  have  k^<<k^  as  t-**”,  and 


k  ^  ^  m  k  as  t 
w 
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As  t-*»,  the  Michell  thin-ship  approximation  k^  and  the  Hogner  flat-ship  approxi¬ 
mation  k^j  become 

2  2 

kjl  ^  F  (l-cosa)/dt  as  t 
2  2  3 

kH  ^  2F  (sina) /b  t  as  t  -*• 

We  thus  have  k^<'<k^<<k^  in  the  limit  t-*00.  However,  the  waterline  integral  k^t 

•) 

and  the  Michell  approximation  k^  are  of  the  same  order,  namely  l/t“,  as  t  . 

As  a  matter  of  fact,  the  above-given  large  t  asymptotic  approximation  of  the 
Michell  function  k^  may  be  seen  to  correspond  to  the  thin-ship  limit  b=0  of 
the  large  t  approximation  of  the  waterline  integral  k^.  Finally,  the  Maruo-Tuck- 
Vossers  approximation  k  vie  Ids 

kMTV  1  as  c 

so  that  the  Havelock  wave-resistance  integral  would  in  fact  be  divergent.  Indeed, 
the  derivative  A' (x)  of  the  ship  cross-sectional  area  A(x)  must  be  continuous  and 
vanish  at  the  bow  and  stern  for  the  wave  resistance  to  exist.  This  condition  is 
extremely  restrictive,  and  in  fact  is  not  satisfied  for  usual  ship  forms.  If  the 
beam/length  ratio  b  and  the  draft /length  ratio  d  are  small,  of  the  same  order  of 


magnitude  e  say,  the  functions  k^,  k^°\  k^ 

,  _2 ,  n  n+1 

k  'v  F  he  t  as  t  -*■  », 


and 


are  of  the  form 


(0) 


where  n=l  for  the  waterline  integral  kw,  the  zeroth  approximation  k  and  the 
Michell  approximation  k^,  n=2  for  the  Hogner  flat-ship  approximation  k^ ,  and  n=3 
for  the  Hogner  approximation  k^ . 

In  the  limit  t=0,  corresponding  to  transverse  waves,  the  Maruo-Tuck-Vossers 
approximation  k^y  and  the  Hogner  flat-ship  approximation  k^  are  identical, 
given  by 

k^  =  k^  =  l-2F2sin(l/2F2)  for  t  =  0. 

The  Michell  approximation  k^  and  the  Hogner  approximation  k^  are  also  identical, 
and  given  by 

1^=1^=  F2 [ { 1-1/exp (d/F2) } /d-2s in (1 /2F2)+Ad{ 1-cos (1/2F2) } ] / (l+4d2)  for  t  =  0. 

This  expression  for  k^  and  k^  may  be  seen  to  become  identical  to  the  expression 
for  k^j  and  k^.^  in  the  limit  d>0.  However,  this  limit  is  obviously  not  uniform 
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.  2 
with  respect  to  the  Froude  number.  In  particular,  we  have  k=l  if  d<<F  <■<!, 

2  2  H 

whereas  we  have  k^-F  /d-0  if  F  <<d<<l.  In  the  limit  t=0,  the  waterline  in¬ 
tegral  k^  is  given  by 

=  -4dF2{l-cos(l/2F2) }/(l+4d2/b2+4d2)  for  t  =  0. 

We  thus  have  k^<<k^  for  t=0,  and 
k^^  'v  k  as  t  -*•  0. 

H 

The  fact  that  k^  ■  k^  as  t->oo,  even  though  k^<<k^  for  t=0,  emphasizes  the 
importance  of  the  waterline  integral  in  expression  (7.1)  for  the  zeroth-order 
slender-ship  approximation.  It  may  be  interesting  to  examine  in  passing  the 
effect  of  using  the  thin-ship  approximation  y=0  in  the  exponential  term  E  in 
the  expression  for  the  waterline  integral,  which  then  becomes 

f  exD{-iF-2(l+t2)1/2x}n2t  d l. 

’  c  *  y 

This  "thin-ship  waterline  integral"  is  of  the  same  type  as  the  Maruo-Tuck-Vossers 
integral  (7.6).  We  thus  may  expect  the  Havelock  wave-resistance  integral  to  be 
divergent  in  general  (except  for  hull  forms  having  a  smooth  waterline  with  cusped 
ends)  if  the  thin-ship  approximation  is  used  in  the  waterline  integral.  For 
instance,  for  the  simple  planar  hull  examined  previously,  the  function,  k^  say, 
corresponding  to  the  above-defined  "thin-suip  waterline  integral"  is  given  by 

k^  =  -4b2dF2(l-cosa)/(b2+4d2+4b2d2) , 

for  which  the  Havelock  integral  (6.12)  is  divergent. 

The  above  simple  example  thus  demonstrates  that  the  use  of  the  thin-ship  or/and 
flat-ship  approximations  may  have  important  effects.  In  particular,  these  approx¬ 
imations  are  clearly  not  uniform  with  respect  to  the  Froude  number  and  the  integra¬ 
tion  variable  t  in  the  Havelock  integral  (6.5).  Use  of  such  approximations  evidently 
is  inherent  to  the  methods  of  solution  based  on  a  systematic  thin-  or  flat-ship  per¬ 
turbation  scheme,  due  to  the  necessity  of  transferring  the  boundary  condition  at 
the  hull  surface  to  the  ship  centerplane  or  waterplane,  respectively.  It  was  pre¬ 
viously  discussed  in  Noblesse  (1976)  that  this  transfer  of  the  hull  boundary  condition 
mav  be  overly  and  unnecessarily  restrictive,  and  may  in  fact  be  far  more  limiting 
than  the  linearization  of  the  boundary  condition  at  the  sea  surface.  An  essential 
feature  of  the  present  slender-ship  theory  is  precisely  to  avoid  transferring  the 
hull  boundary  condition,  whicli  is  imposed  on  the  exact  hull  surface. 
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8 .  The  first-order  slender-ship  low-Froude-numher  approximation 

The  next  approximation  in  the  sequence  of  slender-ship  wave-resistance 
approximations  is  the  first-order  slender-ship  approximation,  r^^  ,  obtained 
by  taking  the  potential  $  in  expression  (6.6)  for  the  Kochin  free-wave  amplitude 
function  K(t)  as  the  first-order  slender-ship  potential  given  by  equation 

(5.2a).  This  first-order  slender-ship  approximation  will  be  considered  in  the 
following  section.  In  the  present  section,  we  briefly  examine  the  simpler  approx¬ 
imation  obtained  by  taking  the  potential  <p  as  the  zero-Froude-number  limit 


of  the  potential  4> 
/.  _ .  (  ^ _  f  r  r  /  ?- 


The  potential 


is  given  by 


(8.1) 


4w*<1)(t)  =  -/[{ (C-x)<"+(n-y)24‘(^-z)2}  (£-x)  2+(n-v)  ~+(r>z)  *”}  i^4']nx(x)da(x)  ,  (8. 

h 

as  may  be  seen  from  equation  (5.4a).  In  the  common  case  of  a  single-hull  ship 
with  port  and  starboard  symmetry,  expression  (8.1)  can  readily  be  expressed  as 
an  integral  over  the  positive  (port)  half  h+  of  the  mean  hull  h. 

The  potential  may  be  regarded  as  a  first-order  slender-ship  approxima¬ 

tion  to  the  (exact)  zero-Froude-number  (double-hull)  potential  (J>q,  as  was  noted 
previously  in  section  5  and  is  specifically  indicated  in  equation  (5.6).  The 
low-Froude-number  wave-resistance  approximation,  r^  say,  associated  with  the 
zero-Froude-number  potential  is  essentially  identical  to  the  low-speed 
wave-resistance  formulas  proposed  by  Guevel,  Vaussy,  and  Kobus  (1974),  Baba  (1976), 
Maruo  (1977),  and  Kayo  (1978),  as  will  be  shown  below.  The  wave-resistance  approxi¬ 
mation  associated  with  the  first-order  slender-ship  approximation  (f^to  the 

zero-Froude-number  potential  4>o  thus  corresponds  to  a  first-order  slender-ship 
(1)  ^ 

approximation,  r)  say,  to  the  Guevel-Baba-Maruo-Kayo  low-Froude-number  approxima- 


It  is  interesting  to  compare  the  low-Froude-number  wave-resistance  approximation 

r/,  ,  the  first-order  slender-ship  low-Froude-number  approximation  r and  the  zeroth- 
cr  (0) 

order  slender-ship  approximation  r  ,  associated  with  the  approximations  <P~<Pn, 

(1)  (0)_  0 

<t>=  ,  and  )>=i>Q  =0,  respectively.  This  comparison  is  particularly  simple  in  the 

special  case  when  the  hull  h  is  a  vertical  cylinder  with  elliptic  waterline  c,  and 

the  nonlinear  free-surface  flux  q  in  expression  (6.6)  is  ignored.  Indeed,  the  wave- 

resistance  approximations  an(l  r£F  can  then  be  shown  to  be  proportional 

to  one  another.  Specifically,  we  have 

r (0)  -  r(,F/(l+b)2,  rJF}  =  rfF(l  +2b)  2/ ( 1+b) 4  , 
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where  b  is  the  beam/length  ratio  of  the  elliptical  cylinder.  The  relative  errors 


;(0)  =  (r!F~r(0))/r('F  and 
:(0)  =  2b ( l+b/2) / (1+b) “ , 


We  thus  have  £  v  =.174  and 


:(1)  =  (r£F-r^)/r£F  then  are  given  by 
(1)  =  2b2[l+b(l+b/2)/(l+b) )/(l+b)3. 
l:^3^  =  ,016  for  b=.l,  and  l^^  =  .306  and  e^3^=.055  for 


b=.2.  The  first-order  slender-ship  low-Froude-number  approximat ion  r,,^  mav  then 

<_  r 

be  seen  to  differ  from  the  low-Froude-number  approximation  r^F  by  only  a  few  per¬ 
cent  for  a  thin  ellipse.  More  generally,  the  second-and  higher-order  slender-ship 

(k) 

low-Froude-number  wave-resistance  approximations  r)  ,  k%2,  associated  with  the 

(k)  " 

sequence  of  slender-ship  iterative  approximations  <j>  defined  by  the  recurrence 

relation  (5.6),  can  also  be  shown  to  be  proportional  to  the  appioximation  r»_,  and 

(k) 

rapid  convergence  of  the  sequence  of  approximations  r»  to  r^F  may  be  proved. 

Details  will  be  given  elsewhere. 

It  is  noteworthy  that  the  relative  error,  between  the  wave-resistance 

approximations  rl3^  and  r»_  is  significantly  less  than  the  relative  error  between 

tt  -tt  QS 

the  corresponding  potentials  and  <j>  .  Indeed,  for  the  present  case  of  an  ellip- 

(1)  u  u 

ticaL  cvlinder  we  have  V  =1  /(1+b)  ,  as  is  shown  in  Noblesse  and  Triantafyllou  (1980). 

u  fn  ns  (n 

This  yields  the  relat  ive  error  n'  =  (<)>n-<D  ' )  /  4)n=b  /  (1+b) .  We  thus  have  n  v  =.091  for 
CD  u  U  U  .  .  ci) 

b=.l  and  =.167  for  b=.2,  which  may  be  compared  to  e '  =.016  for  b=.l  and  ev  =.055 

for  b=.2.  This  suggests  that  evaluation  of  the  wave  resistance  by  means  of  the  Have¬ 
lock  and  Kochin  formulas  for  determining  the  energy  contained  in  the  waves  following 
the  ship  may  be  preferable  to  pressure  integration  over  the  hull.  Specifically,  use 


of  the  first-order  approximation 


to  the  potential  yields  a  Havelock-Kochin 


wave-resistance  approximation  comparable  to  the  wave  resistance  which  may  be  obtained 

by  hull-integration  of  the  pressure  determined  from  the  second-order  approximation 

.  (2) 

"0  • 

Calculations  by  Chen  and  Noblesse  (1980)  for  the  case  of  a  vertical  cylinder 
with  a  waterline  in  the  shape  of  an  ogive  have  shown  that  differences  between  the 
approximations  r^F^  and  r^F  are  somewhat  larger  than  for  an  elliptical  cylinder  with 
same  beam/length  ratio,  but  remain  small,  of  the  order  of  a  few  percent.  Furthermore, 
differences  between  the  approximations  r^3^  and  r^F  are  significantly  smaller  for 
a  slender  hull  than  for  a  thin  cylindrical  hull  with  same  waterline.  This  stems 
from  the  fact  that  the  potential  is  a  better  approximation  to  the  potential 

i  for  a  slender  body  than  for  a  compnrabl e  thin  cylinder;  for  instance,  the  potential 

ID 

tv  is  about  13%  smaller  than  the  potential  f>n  for  an  elliptical  cylinder  with  beam/ 

0  ( i  \  u 

length  ratio  equal  to  .15,  but  is  only  about  3%  smaller  than  for  an  ellip¬ 

soidal  liul  l  form  with  tilt*  same  value  of  the  beam/length  ratio  and  with  draft/length 
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ratio  equal  to  .05.  Calculations  reported  in  Chen  and  Noblesse  (1980)  in  fact 

show  that  the  approximations  r ^  and  r^  are  practically  indistinguishable  for 

ellipsoidal  hull  forms  with  beam/length  ratios  equal  to  .1  or  .2. 

The  low-Froude -number  approximation  K^F  to  the  Kochin  free-wave  amplitude 

function  K  is  defined  by  equation  (6.6)  or  the  equivalent  equation  (6.9),  in  which 

the  potential  <J>  is  taken  as  the  zero-Froude-number  potential  An  interesting 

alternative  expression  for  the  low-Froude-number  Kochin  function  K..  will  now  he 

t  F 

obtained.  We  start  from  equation  (6.3).  By  using  the  hull  boundary  condition 

2 

(5.3b),  from  which  follows  the  relation  n  +t  3d>./3E-n  t  ci_/3d=34  /lx  as  was 

x  x  0  z  y  0  0 

shown  in  section  5,  equation  (6.3)  yields 


K£F(t) 


0/3n-<t>03E/ln)da 


|  c  (E3<f>0/ 3x-<| 


3E/3x)t  dE 
0  v 


-  f 


oEq(^_)dxdy 


(8.2) 


where  E:E(x;t,F^)  is  the  exponential  function  defined  by  equation 
the  potential  <J>q  and  the  function  c  satisfy  the  Laplace  equation 

2  2  2  2  1/2 

i>Q^l/r  as  r=(x  +y  +z  )  ->«  and  E  is  exponentially  small  as  z-»— 

the  Green  identity 


(6 . 3a) .  Both 
in  d,  and  we  have 
®.  We  then  have 


Ih  ( ^ 3  b q /  3n-i>03E/3n)da  =  (a  (E3i)0/3z-<i>03E/3z)dxdy  =  F2|a<i)Q32E/3x2dxdy , 

2  2  2 

where  equation  (5.3c)  and  the  relation  3E/3z+F  3  E/3x  =0  on  z=0  were  used.  By 
expressing  the  term  4>Q32E/3x2  in  the  form  3  (  <}>q3E/  3x— E  3  4>q/ 3x)  /  3x+f32f0/3x2 ,  and 
using  the  identity 


/0  3  (fglE /Sx-Ecfi^/'ix) /"<xdxdv  =  (<f>(-)9E / 3x— E 3d>0/3x)dy  , 

we  may  then  obtain 

F  2/h(E3(J)0/3n-4)03E/3n)da  +  {c(E3<J>0/3x-<j>g3E/3x)tyd<E.  =  JoE32<j>Q/3x2dxdy , 

where  the  relation  dyft^dE  on  c  was  used.  Use  of  this  identity  in  equation  (8.2) 
finally  yields 

< 

K£F(t)  =  /  .E[32<}>0/3x2-q((i,0)  ]dxdy  ,  with  E  =  exp[-iF~2  (1+t2) 1^2 (x+ty )  ]  .  (8.3) 

If  the  nonlinear  free-surface  correction  term  q(4>g)  is  neglected  and  the  mean 
sea  surface  a  is  replaced  by  the  entire  plane  z=0,  equation  (8.3)  becomes  identical 
to  the  low-speed  approximation  that  was  first  obtained,  in  a  different  manner,  by 
Guevel,  Vaussy,  and  Kobus  (1974).  It  was  subsequently  noted  by  Kayo  (1978)  that 
integration  in  Guevel 's  formula  had  to  be  restricted  to  the  portion  of  the  plane 
z=0  outside  the  hull,  that  is  the  mean  sea  plane  a,  for  uniqueness.  Expression 
(8.3)  for  the  low-Froude-number  Kochin  function  K^F  was  previously  derived  from 
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expression  (8.2)  by  Maruo  (1977),  in  a  manner  different  from  that  shown  here. 

2  2 

If  the  linear  term  9  4>q/9x  and  the  nonlinear  term  q C 4>q)  are  grouped  (as  is 
shown  below)  expression  (8.3)  becomes  identical  to  the  low-speed  approximation 
obtained  by  Baba  (1976),  by  using  yet  a  different  approach.  By  making  use  of 
equations  (5.3a,c),  expression  (2.2)  for  the  nonlinear  correction  term  q 
can  be  put  in  the  form 

q($0)=3|v^0|2/3x  +  (94>0/3x)V2*0  _  V- j  V*Q  |  2V$Q/2  , 

where  V  is  the  two-d imens ional  differential  operator  (3/9x,3/3y),  We  may  then 
obtain 

92^0/9x2-q(<)>0)  =  3(1-94>0/3x)O^0/9x-|V((>0|2/2)/9x  -  9  (9<j>0/ 3y)  (3<>0/3x- 1  V$Q  |  2/2) /9y  ,  (8.4) 

which  is  identical  to  expression  (11)  in  Baba  (1976) . 

In  summary,  the  Guevel-Baba-Maruo-Kayo  low-Froude-number  wave-resistance 

approximation  r,>„  has  been  shown  to  be  essentially  equivalent  (for  most  practical 

‘(1) 

purposes)  to  the  first-order  slender-ship  low-Froude-number  approximation  r„  . 

(1) 

The  latter  approximation,  associated  with  the  potential  <p  ,  may  be  regarded  as 

u  (1) 

a  particular  case  of  the  first-order  slender-ship  approximation  r  ,  associated 

with  the  potential  specifically,  the  potential  is  the  zero-Froude- 

(1)  U 

number  limit  of  the  first-order  potential  tj>  .  A  major  difference  between  the 

potentials  <J>^  and  resides  in  the  waves  that  are  present  in  the  potential 

(1)  u  (1) 

4  but  are  ignored  in  the  zero-Froude-number  potential  4>q  .  For  the  present 

problem  of  wave  resistance,  it  clearly  is  appropriate  to  retain  the  waves  in  the 

potential  . 


9 .  The  first-order  slender-ship  approximation 

The  first-order  slender-ship  wave-resistance  approximation  r  is  obtained 
by  taking  the  potential  $  in  expression  (6.6)  for  the  Kochin  free-wave  amplitude 
function  K(t)  as  the  first-order  slender-ship  potential  <J>  =  $  given  by  equation 

(5.2a),  or  equation  (5.8)  in  the  case  of  a  single-hull  ship  with  port  and  starboard 
symmetrv.  By  using  expression  (3.3)  for  the  Green  function,  the  potential  if 
may  be  expressed  in  the  form 

*(1)(?;F2)  =  *Q1}(f)  +  +  4.^1)(t;F2)  , 

where  (C)  is  the  zero-Froude-number  potential  given  by  equation  (8.11,  and  the 

potentials  <tej^(f;F2)  and  $J^(£;F2)  are  defined  below.  The  nonose  i  llatory  near- 
N  (1)  W 

field  potential  is  defined  by 


4ttcP^1)  (C  ;F2)  =  (h(N+l/r')nxda  +  F2{c(N-l/r')n2tyd£  , 


(9.1a) 


where  N=N(|,x;F2)  is  the  function  given  by  the  integral  (3.5)  or  the  alternative 

,2221/2 

integral  representations  (3.5a)  or  (3.6a),  and  r'=[  (f.-x)  +(n-y)  +(Q+z)  I  •  The 
wave  potential  in  equation  (9.1)  is  given  by 

4rr*^1)(S;F2)  =  f^da  +  F2fcWn2tyd£  ,  (9.1b) 

where  WrW(£,x;F2)  is  the  wave  function  defined  by  the  integral  (3.4). 

By  using  equation  (3.4)  in  equation  (9.1b),  we  may  express  the  wave  potential 
in  the  form 

w 

n^1)(t;F2)  =  Imf’°exp[F-2(l+T2)1/2{(l+T2)1/2C+i(C+Tn)}]K^0)(T;F2)dT  ,  (9.2) 

where  the  function  K^^(r;F2)  is  defined  as 

K^0)(t;F2)  =  fc  En2tvdf  +  F_2fh  exp{F~2  (1+T2)  z  lEn^a  ,  (9.2a) 


(9.2a) 


in  which  c_  and  h ,  represent  the  portions  of  the  mean  waterline  c  and  hull  h  between 
9  f-’  2 

the  ship  bow  and  the  plane  x=f_,  and  E  E(x,y;T,F  )  is  the  function  given  by 

E(x,v;-,F")  =  exp[-iF  2 ( 1+r “ ) * ^ 2 (x+iv ) 1 . 

(0)  2 

Comparison  of  equations  (9.2a)  and  (7.1)  shows  that  the  function  Kf  (t;F  )  is 
closely  related  to  the  zeroth-order  slender-ship  approximation  K^(t;F2) 
to  the  Kochin  free-wave  amplitude  function. 

The  nonosc  i  llatory  near-field  potential  vanishes  as  F-*0  (and  thus  is 

N  (i) 

negligible  in  comparison  with  the  zero-Froude-number  potential  for  sufficiently 

small  values  of  the  Froude  number);  indeed,  we  have 
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N+l/r'  =  0[F2/(r')2]  and  F2(N-l/r')  =  -2F2/r'  +  0[F4/(r')2]  as  F-+0. 


It  must  however  be  noted  that  these  asymptotic  approximations  are  not  uniformly 
valid  with  respect  to  r'.  As  a  matter  of  fact,  we  have 


N+l/r'  =  2/r '  -  (2/F2)  [l-(z+0  /  (r  '+|  x-?  | )  ]  +  0(r'/F4)  as  r'-O,  and 
F2(N-l/r')  =  -2[l-r3/ (r'+|  x-C| )  ]  +  0(r'/F2)  on  c  as  r '*-0. 


The  integrand  of  the  waterline  integral  in  expression  (9.1a)  may  thus  be  seen  to 

be  a  continuous  function.  The  function  N  can  be  evaluated  numerically  without 

difficulty  by  using  expressions  (3.5)  and  (3.5a),  as  is  discussed  in  Noblesse  (1978). 

Indeed,  the  integrand  of  the  integral  in  expression  (3.5a)  is  a  continuous  function 

for  -l<t<l  and  0<r'<oo;  the  integrand  of  the  integral  in  expression  (3.5)  is  also 

continuous  for  -l<t<l  provided  r'  is  not  zero  (a  logarithmic  singularity  emerges 

as  r'-*0).  Expression  (3.5a)  is  well  suited  for  evaluating  the  function  N  for  small 
2 

values  of  r'/F  ,  whereas  expression  (3.5)  is  better  suited  for  large  or  moderate 
2  2 

values  of  r'/F  (for  r'/F  greater  than  about  2.5).  In  particular,  the  integral 

in  expression  (3.5a)  vanishes  as  r'-H),  whereas  the  integral  in  expression  (3.5) 

vanishes  as  r'-*  00 .  Furthermore,  the  three-terms  ascending  series  derived  from 

expression  (3.5a)  in  Noblesse  (1978)  can  be  used  for  sufficiently  small  values  of 
2 

r/F  ,  and  the  asymptotic  expansion  associated  with  expression  (3.6a)  and  given  in 

2 

Noblesse  (1975)  can  be  used  for  small  values  of  |y-n|/F  and  large  (or  even  moderate) 

2 

values  of  r'/F  .  For  practical  purposes,  an  approximate  expression  for  the  function 
N  may  actually  be  sufficient.  A  fairly-simple  algebraic  approximation  to  the  function 
N,  obtained  by  combining  the  above-mentioned  ascending  series  and  asymptotic  expansion, 
will  be  given  elsewhere. 

By  using  equation  (9.1)  in  equation  (6.6),  we  may  express  the  first-order  slender- 
ship  approximation  to  the  Kochin  free-wave  amplitude  function  in  the  form 


Kd)  _  K(0) 


K 


(1)  (1)  (1) 

0  +  h  +  K' 


where  the  functions  K 

.(0) 


(0) 


K 


W 

(1) 


K 


(1) 


(9.3) 


.(1)  K(l) 


and  K 


(1) 


are  defined  below. 


The  function  K 
by  equation  (7.1),  that  is  we  have 


0  ’  *Sa  ’  "W  '  a 

corresponds  to  the  zeroth-order  approximation  defined  previously 


.(0) 


(t;FZ)  =  /  En^t  At  +  F  2/,exp{F  Z(l+t^>z]En  da. 


-2, 


x  y 


.(1) 


.(1) 


(9.4) 

The  functions  and  iC*'  correspond  to  the  use  of  the  approximations  and 

(1)  U  N  u 

4>=<j>^  ,  respectively,  in  the  terms  in  expression  (6.6)  that  depend  linearly  on  4>; 

we  thus  have 
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Kg1)  +  =  F~2(l+t2)1/2fcE(i^+F2(tx9i(-/S£-nzty3^/8d)/a+t2)1/2]tyd£ 

-  F  4(l+r2)  f^expiF  2 (1+t2) zlEMn^-i  (n^+tn^) /  (1+t 2)^'/2 ]da, 
where  ii,  =  i|)(x;F2)  represents  the  nonoscillatory  near-field  potential  <t>^^+ 

ry  ON 

In  equations  (9.4)  and  (9.5),  the  function  E=E(x,y;t,F  )  is  given  by 


E(x,y;t,F2)  =  exp[-iF  2(l+t2)1/2(x+ty) ] . 


Equation  (9.5),  with  y  taken  as  the  wave  potential  ,  also  defines  the 
function  .  However,  it  may  be  more  convenient  to  express  the  function 

in  the  form 


K^1}(t;F2)  =  (1/ it)  |  k^1}  (t  ;t  ;F2)dt  , 


(9.6) 


which  may  be  obtained  by  using  expression  (9.2)  for  the  potential  1 ^ .  The  function 

(1)  2  W 

kw  (t; t ,F  )  is  given  by 


k^x;(T;t,F":)  =  F~2(l4M:2)1//2fcE[iiH-F^(tx9ilj/9£-n2ty3t|j/9d)/(l4d:  V/Z]tyd£ 

-  F  ^  (1+t2)  Jhexp{F  2(2+T2+t2)c}EtJ;[nz-i(nx+tny)/(l+t2)1^2)da  ,  (9.6a) 

2  2 

where  the  functions  E=E(£,n;t,F  )  and  ^=^(E,n;r,F  )  are  defined  as 

2  -221/2 
E(C,n;t,Fz)  =  exp[-iF  X(l+t  X(£+tn)J  , 

-nC,n;T,F2)  *  Im  exp[iF'2(l+T2)1/2(C+xn)]  K<0)(t;F2)  , 

with  Kj°\t;F2)  given  by  equation  (9.2a).  Finally,  the  function  in  equation 

•s  CT 

(9.3)  is  associated  with  the  integral  over  the  mean  sea  surface  o  in  equation  (6.6), 
and  is  given  by 

K^X)  (t;F2)  =  /CTexp{-iF~2(l+tZ)1/2(x+ty)}q(({><:i))dxdy  ,  (9.7) 

where  the  nonlinear  free-surface  flux  q(<J>^)  is  defined  by  equation  (2.2). 

If  the  mean  waterline  c,  the  mean  hull  h,  and  the  mean  sea  surface  o  in  the 
vicinity  of  c  are  subdivided  into  small  rectilinear  segments  and  planar  triangular 
elements,  the  waterline,  hull,  and  sea-surface  integrals  can  be  expressed  as  sums 
of  analytical  functions,  corresponding  to  analytical  integration  over  the  elementary 
rectilinear  segments  and  triangles.  The  details  of  this  numerical  implementation  of 
the  first-order  slender-ship  approximation  to  practical  hull  forms  will,  however, 
not  be  examined  in  this  study,  which  is  primarily  concerned  with  the  exposition 
of  the  theory. 


2x1/2, 


(9.6a) 
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For  practical  applications,  the  first-order  slender-ship  wave-resistance 
approximation  r^  is  regarded  as  the  main  result  of  the  present  slender-ship 
theory  of  wave  resistance.  However,  second-and  higher-order  approximations  can 
also  be  defined.  Two  second-order  approximations  indeed  are  briefly  examined  in 
the  following  section. 
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10.  The  second-order  slender-ship  approximation 


(2) 


A  second-order  slender-ship  wave-resistance  approximation  r  can  be  defined 
by  taking  the  potential  <fi  in  expression  (6.6)  for  the  Kochin  free-wave  amplitude 


(2) 


given  by 


(10.1) 


An 


function  K(t)  as  the  second-order  slender-ship  potential  <j> 

*(2)(?)  =  /1)(C)  -  L(|;<f>(1))  -  F2/aGq($(1))dxdy  , 

where  L(£;<j>^)  is  the  linear  transform  of  <f>^  defined  by  equation  (5. 2c), 
interesting  alternative  to  the  above  second-order  potential  may  be  obtained  hv 
seeking  a  solution  of  equation  (5.2),  in  which  the  nonlinear  free-surface  correction 
integral  will  be  neglected,  in  the  form  4>(0  =  k(£)i {)(/,)  k(F,)4> (f.)  ,  where 

k(£)  is  assumed  to  be  a  slowly-varying  function. 

The  term  5  <(> (x)-<(>(£ )  =  k(x)i|)(x)-k(t)if(t)  =  kip-k*^  in  equation  (5.2c) 

may  be  expressed  in  the  form  (k$-kA^)-k^($-!flA)+kA($-!)A)  .  Similarly, 

the  terms  3  ;>/■)(.'  ^  3(k i?)/3£  and  3i>/?d  =  3(kij))/3d  may  be  expressed  in  the  forms 

-  (kv)/3£-k.3$/3£+k.  and  3  (k$) /Sd-k^/Sd+k^^/Sd ,  respectively.  The  equation 


4>(£)  =  iD(1)  -  L(t;<t>) 


(10.2) 


obtained  by  neglecting  the  nonlinear  free-surface  correction  term  in  equation  (5.2b), 
can  then  be  expressed  in  the  form 

k*[  J*+L(£  ;$)  ]  =  4>*  +  k*L (t;$)  -  L(t;k$). 

By  multiplying  this  equation  by  and  using  the  relations  k^ij)^  =  A  and  k$  =  4> , 
we  may  obtain 

«>*[</>*+L(f;>i»]  =  -D*  +  (lj)  -  ^L(t;<P),  (10.31 

which  may  be  regarded  as  an  alternative  form  of  equation  (10.2). 

If  the  (Neuraann-Kelvin)  potential  <f>  were  actually  proportional  to  the  first- 
order  slender-ship  potential  $  H  <|>^,  the  term  <j> *$-$*<!>  would  vanish,  and  equation 
(10.3)  would  yield  the  solution 


4>(l)  =  [^(1)(C)]2/t^(1)(l)+L(C;^(1))]. 

Use  of  the  assumption  jL(£;<f>^)  |  <<  |<f>^  (£) 

4(0  -  i>(1)(t)  -  L(t;<t(1)). 


(10.4) 


in  equation  (10.4)  yields 


It  may  then  be  seen  from  equation  (10.  J)  that  the  potential  <j>(£)  defined  by  equation 
(10.4)  is  approximately  equal  to  the  second-order  potential  <f^2^(£),  if  the  nonlinea 
free-surface  correction  term  is  ignored.  The  potential  ^2>  given  by 
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<t(2)(t)  =  [<fr(1)(f)]2/[<J>(1)(C)+  L(|;<^(1))]  -  F2/  Gq(*(1))dxdy 


(10.5) 


thus  provides  an  alternative,  computationally-equivalent ,  second-order  slender- 
ship  approximation  to  the  velocity  potential  <J>. 

(2) 

In  the  zero-Froude-number  limit,  the  second-order  potential  and  the 

^(2)  (2)  <y(2) 

modified  second-order  approximation  <j>  become  the  potentials  and  <jig  , 

respectively,  given  by 


.  (2) 


(O  = 


U)  -  L0(5;4>0  }  » 


(10.6a) 

(10.6b) 


*02)(f)  =  [<t>Q1)(f)]2/[<t>Q1)(t)+L0(|;<(.Q1))]  ,  (10.6b) 

where  is  the  first-order  potential  given  by  equation  (5.4a)  and  Lg(£;<t>g^) 

is  the  linear  transform  of  <f>g^  defined  by  equation  (5.4b).  It  is  proved  in 

Noblesse  and  Triantafyllou  (1980)  that  the  modified  second-order  approximation 

r'j(2) 

^  is  actually  exact  in  the  particular  cases  of  ellipsoidal  hull  forms.  Further¬ 
more,  calculations  for  two-dimensional  flows  about  ogives  for  various  values  of 
the  thickness  ratio,  have  shown  that  the  potential  provides  a  fairly-accurate 

approximation  to  the  exact  potential  <J>n  (even  for  large  values  of  the  thickness 

U  (2) 

ratio),  superior  to  the  straightforward  second-order  potential  $g  (although  this 

approximation  is  quite  satisfactory  for  values  of  the  thickness  ratio  equal  to 

.1  and  .2).  These  results  for  the  zero-Froude-number  case  suggest  that  the  modified 

<\,(2) 

second-order  slender-ship  potential  ij>  given  by  equation  (10.5)  may  provide  a 
fairly-accurate  approximation  to  the  velocity  potential  <(>. 
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Appendix  A.  Interpretation  of  the  waterline  integral 

Let  us  consider  the  closed  body  consisting  of  the  mean  hull  surface  h  closed 
by  a  horizontal  "lid",  £  say,  slightly  submerged  a  depth  6  below  the  mean  sea  plain 
z  =  0,  so  that  the  "lid"  £  becomes  the  "interior  waterplane"  in  the  zero-sub¬ 
mergence  limit  6=0.  The  integral  identity  corresponding  to  equation  (4.9)  for 
the  fully-submerged  body  h+£  is 

0  9  9  9 

<j>*  =  / jGv  'J’dv  -  J  G(d<p/ 3z+F  3  $/3x  )dxdy  +  f  £  [G?<J>/3n-(4>-4,A)  3G/3n]da  ,  (A.l) 

i 

where  d  is  the  mean  flow  domain  outside  the  body  h+£ .  This  identity  may  be  expressed 
in  the  form 

2  2  2  2 

=  f,G7  $dv  -  f  G(3<|>/jz+F  3  tji/3x  )dxdv  +  f.  [G3(|>/3n- ($-<}> . )  3G/  *n]da  +  I  , 
where  the  term  I  is  given  by 

I  =  -/  G(34>/3z+F232f./3x2)dxdy  +  /^[GSiji/Sz-^-^)  3G/3z]dxdy  , 

i 

since  we  have  3)>/3n  3)>/3z  and  3G/3n  =  3G/3z  on  the  horizontal  lid  £. 

By  expressing  the  integrand  of  the  integral  over  the  lid  £  in  the  form 

G(3^/3z+F232f>/3x2)  -  (4>-4>*)  (3G/3z+F232G/3x2)  -  F23  [G3<j>/ 3x- (<{>-<)>*)  3G/3x]/3x  , 

and  using  the  identity 

/^[G^/Sx-U-l’JSG/lxI/lxdxdy  =  -Jc  [G3ij>/ 3x- (<j>-ij>*)  3G/3x]dy  , 
we  may  then  obtain 

2?  2  •  2  2  2 

I  =  -/  G(3$/3z+F  3^<{/3x  )dxdy  +  j^GO^/Sz+F  3  <p/dx  )dxdy 

-J^ ($-$*)  (3G/3z+F232G/3x2)dxdy  +  F2/_.[G34>/3x-(4>-4>*)3G/3x]dy. 

In  the  zero-submergence  limit  6+0,  the  first  two  integrals  cancel  out  one  another 
and  the  third  integral  vanishes,  so  that  the  above  expression  reduces  to  the  water¬ 
line  integral 

I  =  F2/c  [G3fi/3x-(ii>-'f>#)  3G/3x]dy. 

The  identity  (A.  11  for  the  slightly-submerged  closed  body  h+£  thus  becomes  the 
identity  (4.9)  for  the  sea-surface  piercing  hull  h,  and  the  waterline  integral  in 
equation  (4.9)  may  he  seen  to  stem  from  the  effect  of  the  lid  £  closing  the 
slightly-submerged  hull  h. 
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Appendix  B.  Integral  identities  for  the  interior  problem  and  for  the 
combined  exterior-interior  problems 

The  "interior  potential"  p1 ,  defined  in  the  interior  domain  d^,  can  he  shown 
to  satisfy  the  identities 

=  f,  GV2<j,1dv  -  (  G(3<i,1/Sz+F2a2<i>1/Sx2)dxdy  -  f,  (G3<t>1/3n-(i.19G/an)da 

*  d .  o .  h 

i  2  .  i 

-  F  jc(G34>1/3x-<}>13G/3x)dy  (B.la) 

in  the  interior  domain,  that  is  for  Z  in  d.+o.-h-c  , 

ii 

0  =  f.  GV2(J>Ldv  -  f  G(3(()1/3z+F232(ii1/3x2)dxdy  -  f,  (GS^1 /Sn-.).1 3G/ 3n)  da 
i  ->  .1 

-  F  /c(G3<j)1/3x-4)13G/3x)dy  (B.lb) 

in  the  exterior  domain,  that  is  for  %  in  d+o-h-c,  and 

$i/2  =  GV2(j.1dv  -  f  G(34>1/3z+F232$1/3x2)dxdy  -  f,  (G34)1/3n-(t)13G/3n)da 

*  d .  o .  J  n 

1  2  '  1  i 

-  F  J  (G3<T/3x-<r3G/3x)dy  (B.lc) 

if  C  is  exactly  on  the  boundary  surface  h+c  (provided  h+c  is  smooth  at  the  point  £). 
The  integral  identity  corresponding  to  equation  (4.9)  takes  the  form 

0  =  J  GV2(J)1dv  -  /  G(3<j>1/3z+F232<J>1/3x2)dxdy  -  Jh[G3<i>1/3n-(<j>1-<j>*)  3G/3n]da 

-  F  /c[G3<|)1/3x-((|>  -<j.p3G/3x]dy.  (B.2) 

This  integral  identity,  like  equation  (4.9),  is  valid  for  any  point  £,  whether  inside, 
outside,  or  exactly  on  the  surface  h+c,  and  indeed  is  equivalent  to  the  set  of  the 
three  classical  identities  (B.la,b,c). 

If  we  add  the  integral  identities  (4.6)  and  (B.lb),  we  may  obtain  the  relation 

=  !\GV2cf.dv  +  f  GV2i1dv  -  /  G(3<j>/3z+F232<J>/3x2)dxdy  -  f  G(3()>1/3z+F232(t>1/3x2)dxdv 
i  i 

+  <’h[G(34>/3n-3<ti1/3n)-(4>-(ti1)3G/3n]da  +  F2/c  [  G  (3  d>/ 3x- 3  <*> X  /  9x)  -  ( 1 )  3G/3x]dy . 

Addition  of  the  integral  identities  (4.6')  and  (B.la)  yields  the  same  relation, 
except  for  the  fact  that  $  on  the  left  side  is  replaced  by  <J>^.  We  then  have  the 
relation 

9  9  9  9 

p.  =  f  _GV  <{) dv  -  /  _G(3<|>/3z+F  3  i}>/3x  )dxdy 

*  '  Z<\J  z=u 

+  [,  [G(3(j>/9n-3(J>1  /  Dn )  —  ( <t>— *  )  <>C/  3n  l<!a  +  F2f  [G(34>/3x-3<|> 1  /.lx)- (<f>-<t> 5  )  3G/ 3x  ]dy ,  (B.l) 

7  h  c 

where  on  the  left  side  and  in  the  first  two  integrals  on  the  right  side  corresponds 
to  i>  or  cj> 1  for  points  outside  or  inside  the  hull  surface  h+c,  respectively.  Naturally, 
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the  integral  relation  (B.3)  can  also  be  obtained  by  adding  identities  (4.9) 
and  (B.2).  Indeed,  this  yields 

(l-C^H*  +  =  /2< gGV  <j)dv  -  Jz_QG(9<j)/8z+F^3^<j)/3x^)dxdy 

+  /h[G(34>/3n-3!))1/3n)-((}>-^1)3G/3n3da  +  F2/c  [G(3<J./3x-3<(>1/3x)-($-<j)i)3G/3x]dy  , 

where  0^=  j ^3G/3nda+F2/c3G/3xdy.  It  was  shown  below  equation  (4.9)  that  we  hav> 
C  =0  or  1  for  f  outside  or  inside  h+c,  so  that  the  expression  (1-C .  )<)>  +C  .4' *  is 

i  ->  li 

identical  to  or  for  £  outside  or  inside  h+c,  and  the  above  relation  is 
identical  to  relation  (B.3). 
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Alternative  expressions  for  the  Kochin  function 


It  was  shown  in  section  8  that  expression  (6.6)  for  the  low-Froude-number 
Kochin  free-wave  amplitude  function  K^  ,  obtained  by  taking  the  potential  <f>  as 
the  zero-Froude-number  potential  ^  ,  could  be  expressed  in  the  Guevel-Baba-Maruo- 
Kayo  form  shown  in  equation  (8.3).  An  expression  for  the  Kochin  function  similar 
to  expression  (8.3)  can  also  be  obtained  if  the  potential  <{>  in  expression  (6.3) 
is  approximated  by  the  first-order  slender-ship  zero-Froude-number  potential 


or  more  generally  by  any  waveless  potential,  ip  say,  that  vanishes  sufficiently 

2  2  1/2 

rapidly  as  (x  +y  )  -*■  to  permit  the  use  of  Green's  identity  for  the  potential 


[p  and  the  function  E  in  the  unbounded  exterior  domain  d. 

2 

By  writing  the  expression  n^+t^Sik/S-C-n^t^iJj/Sd  in  the  waterline  integral 

in  equation  (6.3)  in  the  form  n  (n  -9iJj/3n)  +  n  3i(//3n+t  dip/d£-n  t  dip/dd  ?.  n  (n  -3iJ//3n) 

xx  x  x  z  y  xx 

+  3i|'/3x  ,  we  may  obtain 

K(t)  =  F  2/h(E3W  3n-ii)3E/3n)  da  +  fc  (Edi>/dx-tpdE/dx)t^d-E  -  {^Eq  (ijOdxdy 

+  F  (,  E( n  -3tJ//3n)da  +  f  E(n  -3iJj/3n)n  t  &£.  (C.l) 

;h  x  1  c  x  x  y 

By  applying  Green's  identity  to  the  functions  E  and  \|>  in  the  mean  flow  domain  d,  as 
was  done  for  deriving  equation  (8.3)  from  equation  (8.2),  we  can  finally  obtain 

K(t)  =  F"2{  E(n  -3^/3n)da  +  J  E  (n  -3<J>/3n)n  t  dE  +  F ~2J  E[ 3i^/3z+F232^/3x2-F2q(^)  ]dxdy.  (C.2| 
n  x  c  x  x  y  (J 

If  the  potential  ij>  is  taken  as  the  zero-Froude-number  potential  4>q  ,  we  have  3^/3n=nx 

on  h+c  and  3i|^/3z=0  on  a,  and  equation  (C.2)  becomes  identical  to  equation  (8.3)  for 

the  low-Froude-number  approximation  Kp  .  On  the  other  hand,  if  \p  is  taken  as  the 
(1)  tF 

(  n  1  A  '  *  rtf-  -J  /in  /  P  9  \  »»-{  Jn  hVin  f  1  rtr.Ti  n  r*  o  1  h  Ck  n  f"  -f  1TO  nvnv'/X'  ti  1  rvn  f  A  v  f-ko 


potential 


,  equation  (C.2)  yields  the  following  alternative  expression  for  the 


first-order  slender-ship  low-Froude-number  approximation  K^  : 

K^}(t)  =  F~2/hE(nx-3^i;>/3n)da  +  \E  (n^-d^  /  3n)nxtydE  +  fgE[  32^1)/3x2-q(^1))  ]dxdy. 

Another  expression  for  the  Kochin  function  K(t)  may  be  obtained  if  the  potential 
in  equation  (6.3)  is  approximated  by  the  first-order  potential  or  its  zero- 

Froude-number  limit  »  which  are  defined  in  the  domain  d  inside  the  hull  surface 

U  (l)1  (1) 

h  as  well  as  in  the  exterior  domain  d.  The  potential  4'=4>  or  4>q  is  continuous 

across  the  hull  surface  h+c.  Equation  (C.l)  may  thus  be  expressed  in  the  form 

K(t)  =  F-2/h  (E3ijj/3n-ii/3E/3n)da  +  /c  (E 3iJj/3x-tf/3E/3x) t  di  -  /^Eq(i^)dxdy 

hi  Ci  y 

+  F  J,  E(n  -3i|j/3n)da  +  (  E(n  -3i/j/3n)n  t  d-f  ,  (C.4) 

n .  x  c  ,  x  xy 

i  i 

where  h.  and  c.  represent  the  interior  sides  of  the  mean  hull  surface  h  and  waterline 
c.  In  other  words,  the  derivatives  3ij;/3n  and  3ip/3x,  which  are  discontinuous  across 
h+c,  are  evaluated  on  the  interior  side  of  h+c  in  equation  (C.4),  whereas 
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they  are  evaluated  on  the  exterior  side  of  h+c  in  equation  (C.l). 

Both  the  potential  tp  and  the  exponential  function  E  satisfy  the  Laplace 
equation  in  the  interior  domain  d^  ,  so  that  we  may  use  the  Green  identity 

/,  (Edip/ Bn-ipdE/ 3n)da  =  f  (tl>dE/3z-ES\p/dz)dxdy . 
n .  o . 

1  1  2  2  2 

By  using  the  fact  that  the  function  E  satisfies  the  relation  5E/3 z+F  3  E/3x  =0 

2 

on  a.  ,  we  may  express  the  integrand  i|>3E/3z-E3ij>/3z  in  the  form  F  3  (E3^/ 3x-^3E /3x)/ 3x 

1  2  2  2 

-E( dip/dz+F  3  p/dx  ).  Use  of  the  relation 
f  3  (Edip/dx-ipdE/dx)  / Sxdxdy  =  -/  (ESiJ'/Sx-iJ^E/SxH  d £  , 

J  0  C .  V 

i  i 

then  yields 

2  2  2  2 
/  ( E3’J//3n-i()3E/3n)da  =  -F  f  ( E3ip/3x-iji3  E/3x)  t  dl  -  E(3<Ji/3z+F  3  ip  /  3x  )dxdy. 

hi  Ci  V  °i 

Equation  (C.4)  may  finally  be  expressed  in  the  form 


K(t)  =  F  2J  E(n  — 3ip/ 3n) da  +  /  E(n  -3i|i/3n)n  t  dZ 
n .  x  x  x  y 

-  F~“ f  E(3\J//3z+F232i|>/3x2)dxdy  -  /  Eq(i</)dxdy.  (C.5) 

i 

If  the  potential  in  equation  (C.5)  is  taken  as  the  first-order  slender-ship 
(13  22  2 

potential  4>  ,  we  have  3iJi/3z+F  3  ip/dx  5  0  on  a,  ,  and  equation  (C.5)  becomes 

K(1)(t)  =  F-2}  E(n  -3<j>(L)/3n)da  +  J  E(n  -3*(1)/3n)n  t  d Z  -  /  Eq4(1))dxdy .  (C.6) 

'  n .  x  x  x  y  o 

Equation  (C.6)  provides  an  alternative  form  of  the  expression  for  the  first-order 

slender-ship  approximation  defined  by  equations  (9.3)  through  (9.7).  If  the 

potential  <|>  is  taken  as  the  zero-Froude-number  potential  »  equation  (C.5)  yields 

KEF)(C)  =  F'^h  E(nx-3*o1>/3n)da  +  /.  E{  /3n)Vyd£  -  fg  Ed2^1} /dx2dxdy 

-  /a  Eq  <  < i>Ql^)dxdy.  (C.  7) 

An  alternative  form  of  expression  (C.6)  may  be  obtained  by  using  the  relation 
[3i(1)/3n]h  -  [3$(1)/3n]h  =  nx  , 

where  the  notation  [  ],  and  [  ].  implies  that  the  term  within  the  brackets, 

(1)  n  ni 

namely  3<t>  /3n,  is  evaluated  on  the  exterior  and  interior  sides  of  the  mean  hull 

surface,  respectively.  This  then  yields 


nx  -  [3*(1)/3n]h 


nx  +  nx  -  [3<T  '/9n]h  , 


and  equation  (C.6)  may  be  expressed  in  the  form 

w  (^)  \  _  1/  (0)  /..x  I  /  r  /  _i_  f  r  / 


KU;(t)  =  KW(t)  +  F  /  E(n  -34>U;/3n)da  +  /  E(n  -3$U;/3n)n  t  di  -  /  Eq  (d>  U')dxdy,  (CA 

fi  x  c  x  x  y  g 
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where  equation  (7.1)  was  used.  The  following  alternative  form  of  expression  (C.7): 

K£F)(t)  =  K(°)(t)  +  F  2/hE(V34>01)/3n)da  +  /cE(nx-34'01)/3n)nxtyd'£  ~  fa  2|)>q1)/3x2 
_  /aEq('f>Q1^)dxdy  1 

may  be  obtained  likewise. 
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